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Synopsis 


An important problem in control engineering is to design a single feedback control law which 
will stabilize finite number of plants. This is known as the simultaneous stabilization (SS) 
problem.* This problem is motivated by several practical considerations, for instance, the 
dynamics of an aircraft varies greatly with its altitude and speed and requires stabilization 
at all flight conditions. An industrial plant is required to remain stable m different modes 
due to possible failures of components such as sensors and actuators. Often a linearized 
model of a non-linear system operating at different operating points has time varying 
parameters. In such situations, it is important to design a feedback controller which 
will stabilize the system in spite of such changes. This thesis is aimed at studying the 
simultaneous quadratic stabilization problem. 

Consider a finite collection of Linear Time Invariant (LTI) systems denoted by {A, B l ), 
i = 1,2,. . , l and described by the state equations 

x = A % x + B z u 

where x € R n is the state vector and u G R m denotes control input. The basic problem 
considered in this thesis is to find a single state feedback control law u = — Kx such 
that there is a quadratic Lyapunov function V(x) — x*Px ( P being a positive definite 
matrix) which has negative time derivative along the solution of each of the closed loop 
systems. This problem is referred to as the Simultaneous Quadratic (SQ) stabilization 
problem. More specifically, SQ stabilization problem involves the design of a single control 
law u = — Kx and a positive definite matrix P such that the derivative of a common 
quadratic Lyapunov function V(x) = x*Px will be negative definite for the solutions of 
each of the closed loop system. For a general finite set of systems {A,#*}, the problem 
of showing the existence of K and P is quite difficult and open. In the present thesis, we 
attempt to identify different structures of state space matrices {A T } and {A} for which 



SQ stabilization problem is solvable. The thesis is outlined as follows. 

Chapter 1 discusses the importance of simultaneous stabilization problem in control 
theory with a brief literature review of the methods of controller design for uncertain 
systems. 

In chapter 2, definitions and preliminary results for SQ stabilization are introduced. It 
is then established that if a finite set of systems {A,, B}, with fixed input matrix B , are si- 
multaneously transformable into Controllable Companion (CC) form or more generally 
into the Hessenberg form (called as Hessenberg family) then they are SQ stabilizable. Al- 
gorithms are also proposed to achieve stabilization using a single linear static state feedback 
controller It turns out that the family of systems, which are simultaneously transformable 
into CC form, satisfy the matching condition. However, the Hessenberg family is shown 
to be larger than the family satisfying the matching condition. At the end of this chapter, 
conditions for SQ stabilization of a set of systems {A^B*}, with uncertain input matrix, 
are discussed. 

Chapter 3 reports three new families of systems which are SQ stabilizable. These 
families of systems are called (1) partially commutative (2) partially normal and (3) 
a more general family Each of these families is defined by imposing certain assumptions 
on the state space matrices {A,} and {B;} These families of systems are shown to be 
different from the matched uncertain systems. For each case the existence of a simultaneous 
quadratic stabilizing controller is ensured. 

For state feedback controller, it is necessary to measure all the state variable for feed- 
back. However in practice, the design of output feedback controller is a realistic approach. 
Chapter 4 is devoted to the design of static output feedback controller for the simultane- 
ous quadratic stabilization problem. For this purpose, we consider a finite collection of 
minimum phase systems {A t ,B t ,C} described by 

x = A t x + B z u , y = Cx 

where y is the output vector. These systems are, in addition, assumed to be square systems 
(equal number of inputs and outputs) . 

Two new classes of systems are then defined which are SQ stabilizable by output feed- 
back controller. An algorithm for computing a single static output feedback stabilizing 
controller is also presented. 

Chapter 5 is concerned with the design of a simultaneous stabilizing controller for 


IV 



time delay systems and stabilization of time varying systems. Effectively, the results of 
chapter 2 and 3 are extended for special classes of time delay systems and the existence 
of a single stabilizing controller is proved. More specifically, different class of time delay 
systems are shown to be SQ stabilizable by a single memoryless state feedback controller. 
The stabilization of time delay systems is addressed here using a quadratic Lyapunov 
functional. This chapter concludes with the design of a single stabilizing controller for 
special classes of time varying systems 

Finally, m chapter 6. we take two practical examples for which the solution of SQ 
stabilization problem is applied. These examples are (1) stabilizing controller design for 
an uncertain aircraft model (2) Tracking controller design for a robot manipulator 

In this thesis, parametric conditions (structures of system matrix A t and B t ) are derived 
for which SQ stabilization problem is solvable Chapter 7 summarizes all such parametric 
conditions and the overall conclusions of the thesis. 


v 



Acknowledgements 


I wish to express my deep sense of gratitude to my supervisor Dr. V. R. Sule for his 
guidance, advice and encouragement. It was a rewarding experience for me to interact 
with him and gam insight into the topics of research. The various values that I tried to 
learn from him shall remain a source of inspiration for me forever. I owe it to my teachers 
particularly professors A. Ghosh, S S Prabhu, V.P. Sinha, R.N. Biswas, P.R.K. Rao, B. 
Sarkar, P. Sircar, K.E. Hole, A.K Rama for the deep insights given through the various 
courses they taught. 

During the course of this thesis, a significant fraction of time has been spent m the 
company of my various friends. We have discussed a lot of events, issues, ideas and people 
threadbare. Together we have utilized and wasted quite a lot of time. I have greatly 
benefited from the cooperation and help extended by them. In this regard I am thankful 
to Joy, Pavitra, Arun, Sandeep, Nilu, Susa, Anirban, Biswajit, Joydeep, Deepak, Prem- 

4 * 

ananda, Soumya and Shymal. I am also specially thankful to Anirban and Joy for their 
willingness to discuss various technical ideas with me. 

I have benefited greatly from the interaction with many research scholars. For this, I 
would like to thank Shafi, Deshpande, Shivnarayan, Chaitanya, Biswarup, Kishore, Pu- 
ranjoy, Sudipta, Balbindar, Babu, GK, Dandapat, Ramesh, Vinod, Jitendra, Tsr, Apu, 
Harish, Rahul. 

I am gratefull to Metia da and Pal da for their homely affection and enquires about me. 

It was impossible to venture into the Ph.D programme without the support and en- 
couragement of my family members. They have willingly taken up difficult tasks to make 
me feel comfortable. I thank them for patiently bearing my repeated assurance that thesis 
will be completed in notime. 


Indra Narayan Kar 


vi 



DEDICATED to 
my Parents 



Contents 


Synopsis iii 

Acknowledgements vi 

List of figures xi 

List of symbols and abbreviations xii 

1 Introduction 1 

1.1 Purpose of the Thesis . ... . . . . . . 2 

1.2 Simultaneous Stabilizing Controller Design 2 

1.2 1 Previous Approaches to Simultaneous Stabilization Problem .... 3 

1.2.2 Quadratic Stabilization 3 

1.2.3 Other Approaches to the Stabilization of Uncertain Systems .... 6 

1.3 Motivation and Problem Considered in this Thesis 7 

1.3.1 Organization of the Thesis 9 

2 Simultaneous Quadratic Stabilization Problem 11 

2.1 Matrix Inequalities Derived in [6] 11 

2.2 Problem Formulation 12 

2.2.1 Discussion of the Definitions and Implications 13 

2.2.2 Coordinate Dependence of SQ Stabilization Problem 14 

2.2.3 Problem Statement 16 

2.3 Preliminaries 16 

2.4 Basic Inequalities of SQ Stabilization Problem 19 

2.5 Stabilization of Simple Systems 23 

viii 



2.6 Stabilization of Matched Systems 23 

2.6.1 The class of simultaneously transformable systems 26 

2.7 Stabilization of Hessenberg Family 27 

2.8 Results For General Input Matrix Case 31 

2.9 Controller Design Algorithm 33 

2.10 Conclusion 35 

3 Stabilization of Partially Commutative and Normal Systems 37 

3 1 SQ Stabilization of Commutative and Normal Systems 37 

3.2 SQ Stabilization for More General Class of Systems 43 

3.3 Stabilization of a Continuously Parameterized Family of Systems 45 

3.4 Conclusion 46 

4 Simultaneous Stabilization by Output Feedback 47 

4.1 Introduction . . 47 

4.2 Problem Definition and Results . . ... 48 

4.3 Stabilization of Matched Systems by Output Feedback . . . 54 

4.4 Results for Mismatched Systems ... . 56 

4.5 Controller Design Algorithm . 57 

4.6 Conclusion . 59 

5 Stabilization of Time Delay and Time Varying Systems 60 

5.1 Preliminaries 60 

5.2 Stability of Time Delay Systems 62 

5.3 Main Results 65 

5.4 Controller Design for Time Varying System 68 

5.5 Conclusion 70 

6 Applications to Case Studies 71 

6.1 Robust Stabilization of an Aircraft 71 

6.2 Controller Design For Robot Manipulator 74 

6.3 Preliminaries 74 

6.4 .Application to Robot Control 75 

6.4.1 Robot Dynamics • Linearized Model 76 


IX 



6.4.2 Controller Design 77 

6.4.3 Numerical Example 80 

6.5 Conclusion 81 

7 Conclusions 84 

7.1 Conclusions 84 

7.2 Scope of Further Research 86 

A Long Proof 87 

A.l Proof of Lemma 2.3 .' 87 

A. 2 Proof of Theorem 2.5 88 

B Long Proof 91 

B. l Proof of Theorem 5.1 91 

C Quadratic Form of a Matrix 93 

Bibliography 94 


x 



List of Figures 


1 Distribution of eigenvalues of A(r) — HK with variation in r 73 

2 Structure of Robot Control Systems 78 

3 Response of first joint . . 82 

4 Response of second joint .... ... 82 

5 Joint angle errors 83 



List of Symbols and Abbreviations 


0 ■ The zero matrix or vector of appropriate dimension 
A < 0 : A is a negative definite matrix. 

A > 0 : A is a positive definite matrix. 

A > ( <)B : A — B is a positive (negative) definite matrix. 

I m : (m x m) identity matrix. 

W : A positive definite matrix. 

A m [A\ : Maximum eigenvalue of the matrix A. 

A m [A] : Minimum eigenvalue of the matrix A. 

A [A] . Eigenvalue of the matrix A. 
max{7i} : Maximum of the set of scalars r y l . 

A~ l ■ Inverse of matrix A. 
det(A) : Determinant of matrix A. 

A* : Complex conjugate transpose of matrix A. 

Ill'll • The norm of matrix E i e., ||i?|| = A h[E*E]. 

V % : for all i. 

{Ai, B t } ■ A set of systems. 

(A t , B t ) : A single system. 

C~ : Open left half of the complex plane. 

R n : n dimensional Eucledian space. 

Re a : Real part of the complex number a. 

SQ : Simultaneous Quadratic (or Simultaneously Quadratically) . 
QS : Quadratic Stabilization. 

SS : Simultaneous Stabilization. 

CC . Controllable Companion. 

LTI : Linear Time Invariant. 

MIMO : Multi-input Multi-output. 


xii 



Chapter 1 
Introduction 


The feedback mechanism has proved to be a powerful method in controlling the response of 
a system m a desired manner. One primary reason for using feedback control is to reduce 
the sensitivity of a system due to model uncertainties and external disturbances which can 
not be modeled exactly. The higher the sensitivity of the system, the larger the expected 
deviation m its output, which is undesirable. The search for ways to design a system 
insensitive to parameter variations, in fact, led to the discovery of electronic feedback 
amplifiers. In a stable closed loop system, the sensitivity is bounded (but not necessarily 
minimum). Hence one of the important qualitative problems in feedback control theory is 
to design a stable closed loop system despite the model uncertainties and perturbations. 

In classical feedback theory, gain and phase margin were used as indices of the stability 
robustness against model perturbations. In the modern work on feedback stabilization 
theory, new problems were posed to consider explicitly the uncertainties and perturbations 
of the system being stabilized. These problems are popularly known as 

1. Robust stabilization under unstructured uncertainties 

2. Robust stabilization under structured (parametric) uncertainties and 

3. Simultaneous stabilization 

Unstructured uncertainty is specified as a norm bound of the model transfer function 
while parametric uncertainty is specified by a family of model transfer function whose 
structure is fixed but the parameters take values over specified intervals. The uncertainty 
of simultaneous stabilization is specified by a finite set of systems, and the problem is to 
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design a controller such that the closed loop system is stable, no matter which system in 
this finite set represents the plant. 


1.1 Purpose of the Thesis 

The purpose of this thesis is to study a special variety of the linear simultaneous stablization 
problem called the Simultaneous Quadratic (SQ) stabilization problem. This thesis is 
aimed at determining the parametric structures of finite number of state space models of 
control systems which can be stabilized by a single controller. We show the usefulness of 
these results for the stabilization of special classes of time varying systems as well as for 
the stabilization of uncertain time delay systems. We include case studies as applications 
of these results to (1) Controller design for robot manipulator for tracking purpose and (2) 
non-switching controller design for uncertain aircraft model. 

1.2 Simultaneous Stabilizing Controller Design 

To motivate the problems of this thesis, we consider a brief review of the simultaneous 
stabilizing controller design problem in control theory. Simultaneous stabilization problem 
is important to many practical situations in control theory. In practice a plant undergoes 
changes in several ways. For instance, a plant with a nonlinear model may have multiple 
operating points. The linearized dynamics of the plant at these points is thus different. Dif- 
ferent operating points also arise due to hardware and software failures, most importantly 
sensor failures. In all the above cases, fixed stabilizing controller for different operating 
points is a desirable choice. 

Another important application of simultaneous stabilization problem is to tackle system 
uncertainty due to uncertain physical parameters or limited information about the plant 
model. The uncertain physical parameters may consist of masses, inertias, spring constants, 
reaction rates, aerodynamic coefficients etc., that are required in mathematical descriptions 
of the plant. Linearization, model order reduction and neglected coupling terms are some 
of the examples of model inaccuracy Any model, therefore, is at best an approximation 
of reality Of course, the specific description of the source of uncertainties depends on 
the physics and the engineering of the particular system in question. Thus an important 
requirement to be satisfied by a controller is the invariance of the closed loop stability 
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property under model uncertainties. 

In view of above considerations, the simultaneous stabilization problem is relevant in 
the following control problems. 

• Design of a fixed stabilizing controller for uncertain system. 

• Controller design with integrity, i.e., controller continues to operate in the presence 
of hardware and software failures. 

• Single stabilizing controller design of 

— An aircraft for varying altitude and speed. 

— A power system for varying load conditions. 

— A missile for varying inertia and altitude. 

1.2.1 Previous Approaches to Simultaneous Stabilization Prob- 
lem 

The simultaneous stabilization (SS) problem has a long history Several important re- 
sults are reported on this problem in [57, 71, 23, 24] using the systematic approaches of 
factorization theory [72]. In [71], Vidyasagar and Viswanatham derived the existence of 
simultaneous stabilizing controller for the case of two plants. Roughly speaking, two plants 
are simultaneously stabilizable iff their difference plant is stabilizable by a stable controller. 
It was also shown that simultaneously stabilizing (N + 1) plants is equivalent to the simul- 
taneous stabilization of N plants with a stable controller. In [81], Youla et al. studied the 
stabilization of a single plant by a stable controller using the property popularly known as 
parity interlacing property (p.i.p) A plant is said to satisfy p i.p. if between any two real 
right half plane zeros there exists even number of right half plane poles. So using p.i.p. 
condition, one can completely solve the simultaneous stabilization problem for two plant 
case. However, the problem of simultaneous stabilization for 3 or more plant cases is still 
open. 

1.2.2 Quadratic Stabilization 

Another approach to simultaneous stabilization problem can be the efficient use of Lya- 
punov stability theory Lyapunov stability theory is conventionally used more for the 
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analysis of linear and nonlinear, time invariant and time varying systems than for control 
system design. In recent times, robust controller design problem for uncertain system is 
addressed using Lyapunov stability theory. Also a concept called Quadratic Stabilization 
(QS) is introduced in modern control literature for studying stabilization of an uncertain 
system represented by the state space description. The basic philosophy of quadratic sta- 
bilization is to ascertain the existence of a single quadratic Lyapunov function to test the 
stability property of the uncertain system. An advantage of this approach is that these 
results can be applied for the stabilization of linear time varying systems in addition to 
the uncertain LTI systems. Consider a linear time varying uncertain system 

x — A(r(t))x -I- B(s(t))u = (A 0 + A A(r(t)))x + ( B 0 4 - A£(r(f)|u (1.1) 

where A A(r(t)), A B(r(t)) represent the system uncertainty and (A 0 , Bo) is the given 
nominal system model. The uncertain parameter vectors r(t) belongs to some prescribed 
bounded set The problem of quadratic stabilization is to find a feedback controller such 
that the closed loop system is stable with a fixed (uncertainty independent) quadratic 
Lyapunov function. In this direction, important results are reported in [42, 4, 35. 54, 
68, 51]. In [4], Barmish gave the necessary and sufficient conditions for the quadratic 
stabilizability and also gave the synthesis procedure for a nonlinear feedback controller. 
But these conditions are not easy to check. 

The Riccati equation approach to the stabilization of uncertain systems has been de- 
veloped by many authors [35, 54]. In this direction, attempt has been made to solve linear 
quadratic stabilization problem where control law is searched m the class of linear feed- 
backs. Along this line, Khargonekar et al. [35] have given some interesting results for the 
QS They consider the following norm bounded time varying uncertainty 

[A A{t) A B(t)) = DF(t)E 

where D.E are given constant matrices and F(t) is such that F*(t)F(t) < /(identity 
matrix) . Then the existence of static state feedback controller for QS problem is associated 
with the existence of a positive definite solution of a parametric Riccati equation. But it 
is not immediately clear for which class of uncertain systems a stabilizing controller will 
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indeed exist. However, there are a number of results which consider the QS of the matched 
uncertain system (the perturbation A A and A B belong to the range space of nominal 
input matrix B 0 ). The uncertain system (1.1) is said to satisfy matching condition if 
the following conditions hold 

AA(r(f)) = B 0 D{r{t)), A ,B(r(t)) = B 0 E(r(t )) 

where ||£'(r(t))|| < 1 for all possible uncertainties r(t). 

Under this assumption, it is known in the literature that the matched uncertain system 
is quadratically stabihzable [32, 49]. Hence, in order to solve QS problem completely, an 
important aspect that remains to be investigated, is to determine whether there exists a 
controller for an uncertain system when the matching condition will not hold. 

K. Wei [77] considered a single input time varying uncertain system where each un- 
certain element varies independently. Then he defined certain geometrical pattern called 
antisymmetric stepwise configuration, by augmenting system matrix A(r(t)) and input ma- 
trix B(r(t)) which were used for the QS of uncertain system by state feedback controller. 
In this case, it is assumed that the system must contain fixed number (equal to the system 
order) of sign invariant entries in proper locations. 

Deviating from the norm bounded uncertainty representation, Geromel et al. [6] consid- 
ered parametric uncertain system defined in convex bounded uncertainty domain. Consider 
a linear uncertain system described by 

x = Ax + Bu (1.2) 

where A and B are matrices belonging to the uncertainty domains D a and Db respectively. 
The uncertainty domains D a and Db are defined as follows : 

N N 

D a = {A: A = fj-i > 0, !} 

1=1 1=1 

Af M 

Db = {B . B = Bji — 0) oq = 1}. 

j=i j=i 

The matrices {A t } and {B,} represent the extreme points of the the convex sets D a and 
Db respectively The collection of all extreme plants {A l ,B J } for i = 1,2,..., A7 and 

j = 1, 2 M corresponds to a set of linear systems. The quadratic stabilization problem 

of uncertain system (1.2) is posed as a joint search problem for a quadratic Lyapunov 
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function and for a single state feedback gain for the collection of all extreme plants. Then 
authors proved that the convex bounded uncertain system is QS by a linear state feedback 
control law iff there exists a positive definite matrix W and a matrix Z 6 R mXn such that 
for all i = 1 , 2 , ... N and j = 1 , 2 , . M, the following matrix inequalities hold 

AW + WA* < B 3 Z + Z*B]. (1.3) 

These matrix inequalities are linear with respect to two unknown matrices (one of which 
must be a positive definite matrix) By exploiting this linearity structures, they used linear 
programming to solve the QS problem. It is noticed that the above matrix inequalities 
are derived in terms of unknown matrices W and Z in addition to the given state space 
parameters {A} and {B l }. However, in practice, it is always advisable to derive the 
existence conditions in terms of given matrices A x and B % 

1.2.3 Other Approaches to the Stabilization of Uncertain Sys- 
tems 

In this section, a few more approaches to the design of robust controller are briefly dis- 
cussed Two general kinds of measure used in the literature to quantify the uncertainty of 
the plant are unstructured and structured uncertainty In the unstructured uncertainty rep- 
resentation, the perturbed plant is represented by the transfer matrix G(s) = G 0 4- A G(s) 
where Go is the nominal transfer function. It is also assumed that for any proper stable 
transfer function r(jw), || AGQubHoo < \r(j w)\ for all frequency w where UAGHcc is the 
maximum singular value of AG. Two types of problem arises m this framework, namely 
analysis and design. In the analysis problem, given a stabilizing controller, one studies the 
extent of uncertainty to which the perturbed system will remain stable. More specifically, 
suppose C(s ) stabilizes the nominal plant, then the perturbed system will remain stable 
iff [18] 

\\C{I + GoC)- 1 r{jw) IV w. 

In the design problem, given a norm bound on the perturbation AGQ'io ), one determines 
whether there exists a robust stabilizing controller. However, this type of norm bounded 
uncertainty with no additional restriction on AG consists of a very large class of perturba- 
tions, and the controller design based on this often leads to a conservative design. In this 
context one may recall that several scientists such as Routh, Hurwitz, Bode. Nyquist and 
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several others, attempted to solve the stabilization problem. In those days the controller 
design was based on the assumed nominal model of the system to be controlled. The main 
focus was to find a stabilizing controller for the nominal plant, with no consideration of 
system uncertainties. 

In many practical processes, the perturbations have a special form because they may 
originate from the variations of physical parameters. In such perturbations, which are 
called structured uncertainty, the transfer matrix or the state space model depends on 
the real uncertain parameter vector p (= po + A p) with nominal value p 0 and the range of 
perturbation A p. The measure of robustness, known as stability robustness margin, is-then 
defined by ||AF|| for which closed loop stability will be preserved. The Routh-Hurwitz 
criterion, the Kharitinov theorem and related results [36, 9] are proved to be valuable 
application tools for the design of robust feedback stabilizing controller. These results are 
used to develop robust stability criteria using the uncertain characteristics polynomial of 
the system. For more details of the existence conditions of stabilizing controller see the 
references [7, 5, 76, 75]. However, m this direction, most of the computational algorithms 
to design controller are developed based on the assumptions that the the characteristic 
polynomial coefficient is an affine function of the controller gain. It is easily verified that 
this indeed is the case for all single input systems [5, 75]. However, for multi input systems, 
this assumption is difficult to satisfy [5]. 

Ackermann [1] used a graphical method to map a stability region (a region in the open 
left half of the complex plane) onto the space of state feedback gains for each plant. Then 
the intersection of all the mapped regions in the gain space gives the desired solution. His 
method is not easily applicable when there are more than two gain parameters Schemitm- 
dorf et al. [59] derived a sufficient condition for the existence of a state feedback controller 
which stabilizes a set of single input systems. They also gave controller design algorithm 
when the sufficient conditions are satisfied. 

1.3 Motivation and Problem Considered in this The- 
sis 

In general, two approaches are broadly used to solve robust stabilization problems, namely, 
the frequency domain approach and the time domain approach. In this thesis the latter 
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approach is followed. In the time domain approach, the notion of Lyapunov stability is 
adopted to solve the simultaneous stabilization problem for a finite number of linear systems 
using a constant feedback controller. In this direction, the results of Geromel et al. [6] give 
fairly general conditions for solving the stabilization of convex bounded uncertain systems. 
However, the above conditions are not constructive, i.e., there is at present no constructive 
algorithm to decide whether the matrix inequalities (1.3) can be solved. In other words, 
the question of existence of the stabilizing state feedback controller is still unresolved for 
the systems (1.2). So there is a need to derive the existence conditions as well as the 
tractable computational algorithm for the SS problem. For this purpose, we first consider 
the following set of systems described by state space systems for i = 1, 2, . . . , l 

x = A z x + B t u , y = Cx (1.4) 

where x is state vector, u is input and y is output vector. Next, we consider the following 
mathematical questions. 

• Given a set of LTI system described by {A,, B z }. what are the conditions (parametric 
structures) on A t and B z for which there exists a single controller gain I\ such that all 
the eigenvalues of each closed loop matrix have negative real part, i.e , X(A t — BJ{) G 
C~ ? This problem corresponds to the design of static state feedback controller for 
simultaneous stabilization problem. 

• Given a set of LTI system described by {A,, B x , C}, what are the structures of A, and 
B l and C for which there exists a single controller gam L such that X(A l —B l LC ) G C~ . 
This problem corresponds to the design of a single static output feedback controller 
for simultaneous stabilization problem. 

To give a partial answer of the above questions, we adopt an approach along the line 
of quadratic stabilization. This approach will be referred as Simultaneous Quadratic (SQ) 
stabilization (the formal definition will be given in chapter 2). The problem of simultaneous 
quadratic stabilization is to find a linear constant feedback gain for a finite number of LTI 
systems such that the eigenvalues of each closed loop matrix will lie m the open left half 
of the complex plane. Furthermore, the stabilizing feedback control law is required to 
ensure that each closed loop system is asymptotically stable with a common quadratic 
Lyapunov function. .In other words, SQ stabilization problem involves the design of a 
single state feedback control law u = —Kx and a single positive definite matrix P such 
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that the derivative of a common quadratic Lyapunov function V ( x ) = x*Px is negative 
definite along the solution of each closed loop 1 system (1.4). Unlike the norm bounded 
time varying uncertainty representation of QS problem, the SQ stabilization problem is 
defined for a finite collection of LTI systems. Furthermore, no concept of nominal system 
is used in SQ stabilization problem. With this motivation, we pose the following questions 

• What are the structures of (those beyond the matched uncertainty case) 

that are SQ stabilizable 2 ( i.e., the stabilizing state feedback controller exists)? 

• What are the conditions for output feedback versions of SQ stabilization? ' 

• How can these results be extended to the SQ stabilization of time delay systems? 

• Whether and how these results would help to find a stabilizing controller for time 
varying situation? 

In this thesis we analyze the SQ stabilization problem and answer the above questions. 

1.3.1 Organization of the Thesis 

The content of the present thesis is divided into different chapters as given below. A brief 
discussion of literature is also included in the introduction part of each chapter. 

In chapter 2, we first define the problem statement of the SQ stabilization of finite 
number of LTI systems {.4 t , B z } by state feedback controller. By assuming fixed B matrix, 
the solution of SQ stabilization problem is derived in terms of the solvability of a set of 
matrix inequalities. The next question considered is as to what the structures of {A} and 
{B z } matrices should be for which SQ stabilization problem is solvable. Towards this goal, 
a set of single input second order systems is considered and the necessary and sufficient 
conditions for the SQ stabilization of these systems are derived. A set of matched uncer- 
tain systems (the formal definition will be given later) are shown to be SQ stabilizable 
As a corollary to this result, it is also shown that if a set of systems are simultaneously 
transformable into controllable companion form, then they also satisfies matching condi- 
tions. Next it is proved that a set of systems which are simultaneously transformable into 

1 After plugging the state feedback control law u = —Kx in (1.4), the closed loop systems are represented 
by x = (.4, - B t K)x 

2 In this thesis, the abbreviation SO is used for the words ‘Simultaneous Quadratic’ or 'Simultaneously 
Quadratically' and this will be understood from the context. 
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Hessenberg form, are also SQ stabilizable. This class of systems is larger than the matched 
uncertain systems. Using the realization {A t ,B l } and under some technical assumptions 
on {B ,} matrices, an algorithm for computing a state feedback controller of SQ stabilizable 
systems is outlined. 

In chapter 3, we introduce two new family of systems and refer them as (l) partially 
commutative and (2) partially normal systems. These family of systems are defined under 
certain assumptions on the state space parameters {A,} and {B,}. The existence of a state 
feedback controller for the SQ stabilization problem are then ensured for these classes of 
systems. Next, a more general class of systems is also identified for which SQ stabilization 
problem is solvable These systems exhibit uncertainties different from those matched 
uncertain systems. 

To implement the state feedback controller, we require the knowledge of all state vari- 
ables. However, m practice it is not always possible to measure all the states. So the design 
feedback control law using measurable outputs is a desirable choice. For this purpose, we 
consider the SQ stabilization problem using static output feedback controller. First, we 
derive the general existence conditions in terms of a set of matrix inequalities. Then dif- 
ferent classes of systems are identified for which SQ stabilization problem is solvable b} 
output feedback controller. These aspects are discussed in chapter 4. 

In chapter 5, the results of chapter 2 and chapter 3 are extended for time delay system; 
and also for time varying systems We begin this chapter by discussing the asymptotn 
stability of linear state delay systems using quadratic Lyapunov functional Then th» 
conditions for simultaneous quadratic stabilizability of a set of time delay systems usinj 
state feedback controller is derived. Different set of delay systems are then identified whicl 
are stabilizable by a single state feedback controller. In the second part of this chapter, w 
discuss the exponential stability of time varying systems. Different class of time varym; 
svstems are reported which are stabilizable by a constant state feedback controller 

We consider the applications of SQ stabilization problem in chapter 6 Two practice 
examples are considered for this purpose. First example shows the design of a singl 
stabilizing controller for an uncertain aircraft model The second example concerns wit 
the design of tracking control law for a two link robot manipulator. Here the specifi 
purpose of the control law is that the joint angle will track the prespecified referenc 

trajectories. 

Chapter 7 concludes the overall thesis. 



Chapter 2 

Simultaneous Quadratic 
Stabilization Problem 


For a convex set of parameterized family of state space model, the quadratic stabilizability 
by linear state feedback control is equivalent to the solvability of a finite number of matrix 
inequalities as given in [6] As stated earlier in chapter 1, a numerical procedure using 
linear programming has been adopted to solve these inequalities [6] The work in this 
chapter is motivated by the lack of existence conditions of the solution of these matrix 
inequalities This chapter is devoted towards obtaining those structures of the state space 
parameters of a finite set of systems for which SQ stabilization problem is solvable (the 
formal definition will be given later). 

The chapter is organized as follows . In the first section 2.1, we briefly state the 
quadratic stabilizability results developed in [6]. Next we define the SQ stabilization prob- 
lem and carry out a preliminary analysis and also develop background results for the sequel 
Section 2.3 is then devoted to the solution of the SQ stabilization problem for special cases 
Finally, the solvability of SQ stabilization problem for special structures of the state spaci 
system parameters is investigated. 


2.1 Matrix Inequalities Derived in [6] 

Consider a uncertain system described by the following state space model 


x = Ax + Bu 


(2.3 
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where matrices A and B belong to the convex sets D a and Db respectively. The uncertainty 
domains D a and D b are defined as follows : 

N N 

D a = {A: A = ^ h tA, Hx > 0, = 1} 

2=1 2=1 

M m 

D b = {B : B = 0 i 3 B 3 ^ a 3 — 1 }‘ 

3 — 1 3 = 1 

The matrices {Aj} and {£h} represent the extreme points of the the convex sets D a and D b 
respectively. Each extreme plant (A„ B 3 ) for % = 1, 2, . . . N and j = 1, 2, . . . , M corresponds 
to a LTI system. Then it is proved in [6] that the convex bounded unqertain system is QS 
by a linear state feedback control law iff there exists a positive definite matrix W and a 
matrix Z G i? mxn such that for all i = 1, 2, ... N and j = 1, 2, . . . M, 

a 1 w + wa; < b 3 z + z*b 

If these inequalities are solvable for matrices Z and W, then one can construct a robust 
stabilizing control law u = -ZW~ l x. It may be noted that these matrix inequalities 
are derived in terms of unknown matrices W and Z in addition to the given state space 
parameters {A,} and { B ,}. However, in practice it is always advisable to derive the 
existence conditions in terms of the given matrices A* and B z alone This observation has 
prompted us to define SQ stabilization problem for a set of LTI systems (not necessarily 
belonging to convex uncertain set). This is discussed below 


2.2 Problem Formulation 

In this section, the simultaneous stabilization by fixed state feedback controller is discussed 
For this purpose, we consider a set of LTI system described by state space equation, 

x l '> = A,s« + B,iiW i = 1,2,. .1 (2.2 

where 6 R n is state vector, G R m is input vector. The matrices A t G R nx ' 
and B t G m x n (m < n) are real matrices respectively This set of systems will b 
denoted by {A,, B l } whereas a single system will be denoted by (A t , B x ). Each pair (A t , B h 
is stabilizable. Based on this state space description, we introduce some definitions c 
simultaneous stabilization by static state feedback controller which are stated next. 
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Definition 2.1 : A set of systems {A t ,B t } is said to be Simultaneously Stabilizable (SS) 
by a state feedback controller if there exists a matrix K such that the the eigen values of 
each closed loop matrix (A t — B l I\) lies in the open left half of the complex plane. 
Definition 2.2 : A single system (A t ,B 1 ) is said to be Quadratically Stabilizable (QS) 
by a state feedback controller if there exists matrices W = W* > 0 and K satisfying the 
following matrix inequality 1 

(A t - B l K)*W~ 1 4- W~ l (A l - B t K) < 0. 

Definition 2.3 : A set of LTI systems {A u B l } is said to be Simultaneously Quadratically 
(SQ) stabilizable by a state feedback controller if there exists matrices W = W* > 0 and 
I\ such that the following matrix inequalities hold 

(At - BJ<)*W ~ 1 + W~ l (A z - BJ<) <OVt. 

2.2.1 Discussion of the Definitions and Implications 

Based on the definition 2.3, the first problem considered in this thesis is stated as follows. 
SQ stabilization problem : Find the different structures of matrices A t and B % such 
that a set of systems {A,, B L ) are SQ stabilizable by state feedback controller. 

Before proceeding further, the following well known result is stated based on the notion 
of Lyapunov stability of homogeneous LTI system x = Ax where x is the state vector. 

Lemma 2.1 The matrix A (associated with a LTI system x = Ax) has all eigen values m 
the open left half of complex plane iff there exists a matrix P = P* > 0 such that either of 
the following inequality holds 
(J) A*P + PA < 0 
or 

(II) AW + W A* < 0 where W = P~\ □ 

In such situation, the derivative of the quadratic Lyapunov function V(x) = x*Px 
will be negative along the solution of x = Ax. The matrix P will be called Lyapunov 
matrix 2 . This lemma will be useful in designing a state and output feedback controller for 
a collection of LTI systems. 

iThe matrix IF -1 is deliberately used to simplify the results derived latter. 

2 The matrix W = P~ l will also be called Lyapunov matrix. 
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In the light of lemma 2.1, it is observed that if a set of systems is SQ stabilizable, 
then the stabilizing controller is given by = —Kx^ and the corresponding quadratic 
Lyapunov function for each system is V(x (l) ) = (:r w )*W -1 rr w . 

It may also be noted that if a single system (A t ,B t ) is stabilizable then it is also 
quadratically stabilizable in the light of lemma 2.1. However, simultaneous stabilization 
of a set of systems (as given definition 2.1) does not necessarily mean that they are also 
SQ stabilizable. But one advantage of SQ stabilization will be that one can give some 
existence condition as well as a numerical algorithm. 

2.2.2 Coordinate Dependence of SQ Stabilization Problem 

Before we progress further towards the solution of SQ stabilization problem, we consider 
an important fact. Note that the SQ stabilization problem as stated above depends on the 
state space representation of the system to be stabilized, and the solution ( i.e., the control 
law u ^ = —Kx^ ) depends on which state space coordinates are used for feedback. In 
other words, if a different state space representation is used for system representation and 
the state feedback is taken with respect to the new state space coordinates, then a non- 
stabilizable family may become stabilizable. To illustrate this fact consider the following 
example. 


Example 2.1 : Consider the following two 2nd order state space models of the systems 
to be stabilized 


i (,) = A,x f,) + Bu<*>, i = l,2 


(2 3) 


where A x 


0 1 
0 0 


i A-2 


0 -1 

0 0 


and Bi = B 2 


0 

1 


Consider the problem of obtaining control laws u® = —Kx^\ % = 1,2 with a single 
feedback gain matrix I\ = [k x k 2 ] such that each closed loop matrix A x - BJ\ for i = 

1,2 is stable, i.e., all the eigen values of (A t - BJ<) are in the open left half plane. 

The characteristic polynomials of the closed loop matrices A x — B X K and A 2 - B 2 K are 
s 2 + k 2 s + ki and s 2 + k 2 s - k x respectively. Now using Routh-Hurwitz criterion, it can be 
concluded that 

(a) A\ — B X K is stable for k x > 0, k 2 > 0 and 

(b) A 2 — B 2 K is stable for k x < 0, k 2 > 0. 
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This is impossible. So there does not exist a single matrix K such that the matrix A l —B l K 
is stable (i.e , all eigenvalues of (A* — B X K) are in the open left half plane) for i — 1,2. 
Hence they are also not SQ stabilizable. 

Now consider two new state variable vectors defined by = T t y W for i — 1,2 


where T x = 


1 0 
0 1 


and T 2 = 


-1 0 

0 1 


Using these state space transformations, we 


get the following two new systems 


V { ' ] = A m yM + 1 = 1,2 


(2.4) 


where A nX = T x 1 A X T X 


0 1 
0 0 


, A n 2 — T 2 1 A 2 T 2 — A n 1 and 


Bnl = T x l B x 


0 

1 


, B n 2 — T 2 1 B 2 — B nX 


Now we see that using a constant matrix /\ = [1 1.732], each system is stablizable by 
control law u ^ = — Ky ^ for x = 1,2. It is also observed that there exists a matrix 
. 0.866 -0.5 . 

W = \ such that the following inequality holds for x = 1, 2 

-0.5 0.866 


(A - Bjcyw- 1 + W~ 1 (A l - BJ<) < 0. 


Thus in the new state space coordinates, {A ni , B n } are SQ stabilizable. The above example 
shows that when state space coordinates of the set of systems {A, A} are given, then the 
state space coordinates are fixed and the control law is to be chosen from the very same 
physical state variables. □ 

In view of the above observations, we shall assume henceforth that the state variables are 
fixed for all the systems and denoted by the same notation i.e., x w = x for alls = 1 , 2 ,..., 1. 
Using the notation u W = «Vt, the system represented by (2.1) is given by 

x = A x x + B x u, i = l,2,...l (2.5) 


The SQ stabilization problem is restated as follows. 

Definition 2.4 : A set of LTI systems {A, B t } defined in (2.5) is said to be simultaneously 
quadratically stabilizable by a state feedback controller u = —Kx if there exists a matrix 
W = IU* >0 such that the following inequalities hold 


(A - BJx'YW- 1 + W~\A X - BJ<) < 0 V z. 
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2.2.3 Problem Statement 

Based on the definition 2.4, the statement of stabilization problem is as follows. 

SQ stabilization problem : Find the different structures of matrices A t and B t such that 
a set of systems {A t , B t } are simultaneously quadratically stabilizable by a state feedback 
controller. 

Observe that SQ stabilization problem involves the design of a single state feedback 
controller u = —Kx such that the derivative of a common quadratic Lyapunov function 
V{x) = x*W~ l x is negative definite for the solution of each closed loop systems. 

2.3 Preliminaries 

In this section, we consider a simplified SQ stabilization problem of as stated 

above, in which input matrix B t are same for each system. However, later in section 2.5, 
this restriction on input matrix will be removed. Now consider a family of LTI systems 
{A,, B}, i = 1, 2, . . ,/ described by 


x = A x x + Bu. (2.6) 

We make the following assumptions about the set of systems {A X ,B}. 


Assumptions 2.1 : The matrix B is full column rank and is of the form 


B = 



where f m is m x m identity matrix. This assumption does not restrict generality, for even 
if B is not of full rank, then by redefining the inputs one can have a B matrix of full rank. 
If the B matrix does not conform to the above structure, then there always exists a single 
similarity transformation such that the above assumption will hold. To further clarify this 
point, consider a set of systems {A X ,B} described in (2.5) where B is not of the above 


structure Let T be the state space transformation such that T l B = 


0 


Irr 


. The new 


system will be given by {A m ,B n } where A m = T 1 A l T and B n = T l B. Now suppose 
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there exists a matrix K and W such that the new systems {A ni ,B n ) are SQ stabilizable, 
i.e., the following inequalities hold 

(A m - B n K)*W~ l + W~ l (A ni - B n K ) < 0 V t. 

After rewriting the above inequalities in terms of original set of systems {A l} B), we obtain 

(A ~ BK 0 )*W - 1 + Wy\A t - BK 0 ) < 0 Vi 

where K 0 = KT~ l and W 0 = TWT*. 

Hence if the set of systems {A nu B n } is SQ stabilizable then {A*,#} is also SQ stabiliz- 
able Similarly, it can be shown that if the set of systems [A t ,B) is SQ stabilizable then 
{A m , B n } is also SQ stabilizable. Thus the above assumption does not lead to any loss of 
generality. 

Assumption 2.2 : For the SQ stabilization problem, we assume that the state feedback 
control law u = —Kx is obtained with K in the form K — (y/2 )B*W~ l where 7 is a 
positive scalar and W is a positive definite matrix. The structure of matrix K is justified 
using the following lemma. 

Finsler’s Lemma : Let M be a given n x n symmetric matrix and let N be an m x n 
matrix such that 

x*Mx < 0 

for all nonzero vectors x € R n such that Nx = 0. Then there exists a constant a > 0 such 
that 

M - crN*N < 0. 

Proof : A proof of this lemma can be found in [35, 33, 21]. □ 

The justification of above assumption 2 2 is as follows. From definition 2.4, a set of 
systems {A,, B) is SQ stabilizable if there exists matrices W = W* > 0 and K such that 
for any vector x G R n the following hold 

- BK)*W~ l + W~\A t - BK)}x < 0 V x E R n and V 1. 

=> x* [. AAV ~ 1 + W~ l A*]x < 0 V 1 and x such that B*x = 0. 

The above condition is a set of matrix inequalities. By Finsler’s lemma for each of the 
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above inequalities, there exists a scalar 7l such that 

x*[A^W~ l + W~ 1 A l - j l W~ 1 BB*W~ 1 ]x < 0 V x G R n (2.7) 

& x*W~ 1 [A t W + WA* - < 0 Vx £ R n 

<=> x*[{A t - BK Z )*W~ X + W~ 1 (A l - BIQ*]x < 0 where K t = ( 7 */2 )B*W~ l . 

The above inequality (2.7) is simplified and rewritten as a matrix inequality 

[A t W + WA* - ~uBB*\ < 0. (2.8) 

In this way, we can derive a set of matrix inequalities of the form 2.8. In other words, 
there exists a set of positive scalars { 7t } such that the following inequalities hold . 

[A t W + W A* - 7l ££*] <0 for i = 1,2,...,/ (2.9) 

Next define a scalar 7 = max{ 7t }. Using this value 7 and noting that the matrix BB* is 
a positive semidefimte matrix, the following inequalities will be satisfied. 

[AtW + WA* -jBB*} < OVi (2.10) 

Furthermore, if (2.10) holds, then m the light of lemma 2.1, V(x) = x*W~ x x is a quadratic 
Lyapunov function for each closed loop system (2.5) with control law given by u = —Kx = 
— ( / j/2)B*W~ 1 x where 7 (> 0) is a scalar. Thus it has been shown that if a set of system 
is SQ stabilizable then one controller gain structure will be K = ('y/2)B*W~ l This com- 
pletes the justification of choosing the controller gain structure. □ 
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The following primary lemma gives a necessary and sufficient conditions for the SQ 
stabilization of a set of systems {A, B}. 

Lemma 2.2 A set of systems {A, B} is SQ stabilizable iff there exists a scalar 7 > 0 and 
a positive definite matrix W > 0 such that 

[AfW + WA* t - iBB*] < 0 Vi (2.11) 

Proof : By definition 2 . 3 , a set of systems {A Z ,B} is SQ stabilizable iff there exists a 
controller gam K and W such that 

(A - BK)*W~ l + W~\A r - BI<) < 0 Vi. 

Now using K = (7/2 )B*W~ 1 , the above inequalities are simplified by pre- and post- 
multiplymg both sides by the symmetric matrix W as 

AffV + W A* - 7 BB* < 0 \/i □ 

Application of the above lemma 2.2 follows m the next section. 


2.4 Basic Inequalities of SQ Stabilization Problem 


In this section, we develop matrix inequalities which facilitate the investigation of SQ 
stabilizable family {A t ,B} Before going to actual discussion, the given matrices A are 
partitioned as follows. 


A 1 A 1 
A __ /1 u ^12 

1 A t A x 
/i 21 /I22 

where A' u € A l 12 € A l 22 e R mxm 


( 2 . 12 ) 


Similarly, the positive definite (n x n) symmetric matrix W appearing in the inequalities 
(2.11) is partitioned as 

" Wx W 2 


w 

Wf w 3 

where \\\ € ^ e K n-mxm ) ^ g R mxm _ 


(2.13) 


The following theorem is stated for the SQ stabilization of {A, B} This result will be 
useful for the design algorithm of state feedback controller. 



2 4 Basic Inequalities of SQ Stabilization Problem 


20 


Theorem 2.1 A finite set of systems {Aj, B }, z = 1,2, . . ,1 is simultaneously quadratically 
stabilizable iff there exists a matrix Wi = Wf > 0 and a matrix F such that the following 
inequalities hold 

(All + A\ 2 F)W l + WfiA\ x + A\ 2 F)* < 0 V i. (2.14) 

Proof (Sufficiency) : r , 

Suppose there exists a matrix W\ = Wf > 0 and a matrix F such that the following 
inequalities holds, 

[(A- n + A' u F)Wi + lU(A' u + A' n F)'] < 0 V ). 

In order to satisfy (2.11), for SQ stabilization, we need to find a positive definite matrix W 
and a scalar 7 > 0 such that $ t , 1 = 1, 2, . . , l will be negative definite where is defined 
by 

^ t = [A t W + WA*-'yBB m ]. (2.15) 

For this purpose we construct a positive definite matrix W in the following way By 
using the matrices F and W x , choose W 2 = FW X Then one can select W 3 such that 

[ Wi Wo 1 

W-i - W*Wf l W 2 > 0 which ensures that W = “is positive definite matrix 

I W.f W 3 

Such a choice of W 3 is always possible. For instance, choose W 3 — oJ m where a > 0 and 
q > \m[W 2 W{ 1 W 2 \ where X M [] is the maximum eigenvalue of matrix [.] 

With the above matrices W and F, the matrix is partitioned as 
<}, = *’ u '“’• 2 where ip\, € 

. W fej 

The different submatrices are defined as follows 

<Pn = (A' u + A\ 2 F)Wi + Wi(A* u + A\ 2 F )* 

4>\ 2 = A‘ U W 2 + A\ 2 W z + WM21Y + W 2 (A\ 2 )*. 

022 = Ai w * + a‘ 2 Ws + w;{A i 21 y + wyw 22 y - 7 / m 

In the above equations F = Wi£ H'T 1 


(2.16) 
(2 17) 
(2.18) 



2.4 Basic Inequalities of SQ Stabilization Problem 


21 


The matrix $ t will be negative definite if [38] 

(-0 011 ^ 0 an d (I I) 022 — (0I2)* (011 ) 1 012 ^ 

It is observed that the matrix 0 l u is negative definite for all i by the hypotheses of 
theorem. Also note that 0 l u and 0 l 12 are independent of 7. So the matrices {022} can be 
made negative definite by choosing a sufficiently large value of 7. Since 0* u and d >\ 2 are 
independent of 7, we can choose 7 such that 

022 — (012) (0n) 1 012 < 0 V i. (2.19) 

The value of 7 is computed as follows. Using the expression of 0 22 from (2.18), the above 
inequalities (2.19) are rewritten as 

1 1n, >s,~ (,«,)* (ft,)' 1 V t (2.20) 

where S, = A‘ 21 W 2 + A' 22 W, + W 2 (A 21 )* + W ? ,(A\ 2 )\ 

Now, for each index 1 , choose a scalar 7 t such that 

7 . > A .„[ S , - ( 0 ‘ 12 )' ( ft,)" 1 fl 2 ] 

where A.v/ [.] is the maximum eigenvalue of matrix [ ]. In this way we choose a set of scalars 
{7i} • 


Let the value of 7 = max{ 7 *}, i = 1,2,.. ,Z. For this value of 7 and W, the matrix 
defined in (2.15), will be negative definite for each % — 1 , 2 , . . . ,Z. The controller gain 
is then given by K = ( 7/2 


Necessary : Suppose, there exists a scalar 7 > 0 and a matrix W > 0 such the following 
matrix inequality holds for each i : 


= [A t W + WA* - 7 BB*] < 0 


(2.211 


Substituting B 


partitioned as 


0 


and partitioning W = 


W x W 2 

w 2 * w 3 


, the matrix are 


0ii 

(« 2 ) 


®i2 

022 J 


By using F = W 2 W X \ the different submatrices are 
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defined as in (2.16-2.18). Since the matrix is negative definite for each i, so matrix 0‘ n 
is also negative definite for each i Hence the proof. □. 

Remarks 2.1 : In theorem 2.1, we prove that a set of system is stabilizable by a 

state feedback controller iff there exists a matrix W x and F such that 

(A l n + A\ 2 F) W x + Wx (Alx + A\ 2 F)* < 0 V z. 

However, to verify this condition, the following steps should be followed 

• Calculate a matrix F such that for each i, the matrix (A l n 4- A\ 2 F) is Hurwitz 3 , i.e., 
all eigen values of (A l u + A\ 2 F) lies in the open left half plane To exist such a matrix 
F, the necessary condition is that each pair (A* u , A\ 2 ) must be stabilizable. 

• Calculate W x = W* > 0 by solving the Lyapunov equations 

(>„ + A\ 2 f) w j + w, (A’ u + a\ 2 f)" < 0 V 1 

So to verify the existence conditions in theorem 2.1, it is necessary that each system pair 
(A\ X) A\ 2 ) must be controllable or stabilizable. 

Next result will be stated for the stabilizability of a system comprising the pair (A l n , A\ 2 ) 
defined in (2.12). 

Lemma 2.3 If a system (A,, B) is controllable (stabilizable) then the pair [A\ l2 A\ 2 ) (con- 
structed from the submatrices of A t ) is also controllable (stabilizable) 

Proof : The details of the proof of this lemma is given in Appendix A. □ 

Based on the theorem 2.1, we derive the existence conditions in terms of the solution 
of reduced order 4 matrix inequalities (2.14) as compared to (2.11) But one needs to know 
the structures of {A t , 5} such that they are SQ stabilizable. The rest of the chapter will 
deal with the characterization such a few classes of systems. Before going to the special 
class of systems, the results of SQ stabilization for a set of second order systems will be 
derived in the next section. 

3 Note that a matrix S is called Hurwitz if all its eigenvalues have negative real parts. The terminology 
arises from the fact that if all eigenvalues of S have negative real parts, then the characteristic polynomial 
of S is a Hurwitz polynomial. 

4 In (2.11). the order of matrix inequalities is n x n whereas in (2.14) the dimension of matrix inequalities 
is n — . m x n — m. 
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2.5 Stabilization of Simple Systems 

In this section, we discuss the stabilization of a set of second order systems. For this pur- 
pose, consider a set of single input second order systems {A^B} where matrix A, and B 
are defined as follows : 


A = 


a 


li 


*21 


a 


12 


x 22 


and B 


0 

1 


Theorem 2.2 A finite set of second order systems {A z ,B},i = 1,2, . 
quadratically stabilizable iff either of the following conditions holds 
(I) a\ 2 is sign invariant entries V i 


l is simultaneously 


or 

(II) a\ x < OVi. 

Proof : Using theorem 2.1, the set of systems {A t , B} will be SQ stabilizable iff for any 
scalar F and W x > 0, the following inequalities hold 

[ (all + a l 12 F) W x + Iff! (a* u + a\ 2 F)* ]<0V l . 

Since W x is a positive scalar, the above inequalities are are equivalently written as 


a l n + a\ 2 F < 0 V i. 


( 2 . 22 ) 


Now the above inequalities (2.22) will hold iff 

(I) a\ 2 is sign invariant entries V i 

or 

(II) a\ x < 0 V 2 . □ 


2.6 Stabilization of Matched Systems 

The concept of matched uncertain system is defined in the robust control literature [32, 49] 
for linear systems with time varying uncertain parameters. In this case, the perturbation of 
system matrix lies m the range space of input matrix 5 as described in chapter 1. However, 
our context, we shall call a finite number of LTI systems ([A*, B } as matched systems if 

■’Suppose an uncertain system is represented by x = (To + Aj4)x + Bqu. Then this system will satisfy 
matching condition if A.4 = B 0 D where D is any matrix of proper dimension 
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the submatrices {A* u } and {A\ 2 } (as defined in (2.12) ) are same for all the systems The 
formal definition is as follows : 


Definition 2.5 : A set of LTI system {A Z ,B} will be called matched system if the sub- 
matrices {A* u } and {A* 12 } are same for all the systems. 


Remarks 2.2 : The above definition has been made by noting the special structure of 

0 1 

matrix B — .In fact for matched systems, the difference between any two matrices 

Bn 

(A, — A,) lies m the range space of B 

Defining this, the SQ stabilization results of such systems can be stated as follows. 

Theorem 2.3 A fimte set of matched systems {A t ,B},i = 1,2, described in (2.5) is 
simultaneously quadratically stabilizable. 

Proof : To prove the above theorem we will consider the special structures of matrices 
.4, and B. Based on the definition 2.5, the matrices A l and B for matched systems are 
defined as follows • 

-4-n Mi 


A, 


Ai Al 

- T -21 ^22 


(2.23) 


and B - 


0 

dm. 


where An E R n 


—mxn—m 


, A u E R n ~ mxm , A 22 € R 


We will show that there exists a positive definite matrix W\ and a matrix F for which 
condition given in theorem 2 1 is satisfied. 

From the inequality (2.14) of theorem 2 1, we need to show that there exists an F and 
Wi such that the following inequality is satisfied 


[ (A u + A U F) W x + W x (A u + A 12 F)* ] < 0. (2.24) 

This is a single Lyapunov inequality. It has positive definite matrix solution Wi iff 
(A u + A 12 F) is Hurwitz matrix. Since (A u , A 12 ) is a stabilizable pair ( since [A Z ,B) is 
stabilizable and using lemma 2.3), there exists non-trivial solutions F for which matrix 
( A u + M 2 F) is Hurwitz By choosing such a matrix F, there exists a matrix W x such that 
(2 24) holds. The computation of a \V l is done by solving the following Lyapunov equation 


[ (A u + A 12 F) W x + (A u + A 12 F)* ] = Q 0 
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where Qo is any (n — m x n — m) negative definite matrix. The existence of matrices Wi 
and F completes the proof of the theorem. □ 

As a special case of the theorem 2.3, it is proved that, if each matrix A t , 2 = 1,2,.. , l is 
in block controllable companion form [62] with respect to the same matrix B then {A t , B} 
are SQ stabilizable. 

Theorem 2.4 A finite set of systems {A t , B},i = 1, 2 , . . . , l in block controllable compan- 
ion form is simultaneously quadratically stabilizable. 

Proof • To prove the above theorem we will consider the special structure of matrices A* 
and B as given below, 

0 I m 0 . . 0 

0 0 I m . . 0 

0 0 0 .. I m 

- otQ —ay —a 2 . — o£ 

where r = n/m and A u € An G R n -mxm j 4^ € Rmxm 

Note that for each system the matrices An and A 12 are same It is also noticed that the 
pair (AmAn) is a controllable pair (being themselves in block companion form). Using 
theorem 2.3, we can conclude that {A^B} are SQ stabilizable. This completes the proof 
of the theorem. □ 





Remarks 2.3 : For single input case (m = 1), the matrix A t of (2.25) will be in the 
familiar controllable companion (CC) form and is given by 

0 1 0 ... 0 

0 0 1 ... 0 

A, = . : 

0 0 0.1 

_ ~ Q i “ a 2 • ■ • ~ a h 

This form will be referred in later stage for comparison purpose. Next we shall state the 
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stabilizability results for a set of systems which are simultaneously transformable into CC 
form 

Corollary 2.1 Consider a set of systems {A t ,B} having the property that there exists a 
nonsingular matrix T for which {T~ l AfT, T~^B} are in controllable companion form, then 
the family {A*, B } is simultaneously quadratically stabilizable. 

Proof : First compute matrix T for any one system. Then transform all given systems 
into CC form using matrix T. The proof directly follows from the result of theorem 2.4 
and as follows. 

Let A ct = T 1 A t T and B c = T~ l B where ( A cl ,B c ) are in CC form. Using the result of 
theorem 2.4, there exists matrices W = W* > 0 and K such that 

£*[ (A« — B c K)*W~ l + W~ l {A cl — B C K) ]x < 0 for any x. 

=> **[ (T -1 A S T - T~ 1 BI\)*W~ 1 +W~ l (T~ l A t T - T~ l BK) ]x < 0. 

=* x*T*[ (A t - BK 0 )*W~ l + IT' 1 (.4, - BK 0 ) }Tx < 0 where W~ l = {T~ l YW^T~ l 
and K 0 = I\T~ l □ 

Remarks 2.4 : Suppose we are given a finite number of systems {A t ,B},i = 1,2,. .,/ 
and they are transformable into CC form using a matrix T. First transform all the systems 
into CC form and then design controller gain A' for the transformed systems. So the actual 
controller gain to stabilize all the original systems {A t ,B} is given by I\ 0 = KT~ l . 

2.6.1 The class of simultaneously transformable systems 

In the last section we have proved that whenever a set of plants {A t ,S},i = l,2,...,Zis 
m controllable companion form there exits a single Lyapunov function for the SQ stabi- 
lizability of the whole set. We show here that the class of single input systems {A t ,B}. 
which are simultaneously transformable into CC form, is equal to the class of systems 
satisfying matching condition. Suppose we are given a controllable pair {Ai,B) and T be 
the transformation matrix such that A cX = T~ l A\T,B c = T~ l B. Construct the matrix 
T as T = [ ei,e 2 , . . e„ ] as given in Wonham [79] using the characteristic polynomial of 
matrix „4i Here the column vector e n is nothing but the vector B itself. After defining T. 
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the following class of systems can be transformed into CC form : 

0 ... 0 ’ 

A a — Ax +T T -1 where d = f dx, . . . , d n } is a arbitrary n dimensional row 

0 0 0 L J 

_dx ... d n 

vector. After simplification, A a can be written as A a = A x + BQ where Q is any arbitrary 
row vector. So the perturbation matrix (A Q — A x ) lies in the range space of input matrix 
B and hence satisfies the matching condition. 

2.7 Stabilization of Hessenberg Family 

In subsection (2.6), it has been shown that the matched family of systems are simultane- 
ously stabilizable. To enlarge this class of system, we consider a family of systems which 
will be referred as Hessenberg family. In this section, we shall discuss the stabilization of a 
set of single input systems which are m Hessenberg family For this purpose, we consider 
a set of n dimensional state space system (A^n), B(n)} described as 

x = A l (n)x 4- B{n)u. (2.27) 

where x & R n and u G R. In the above equation, the matrices {A t (n)} are assumed to be 
in Hessenberg form [41], i.e., 

h\x 1 0 ... 0 

h\x h\ 2 1 ... 0 

A,(n) = ! : ! i € R nXn (2.28) 

Ki-1,1 ^n— 1,2 ^n— 1,3 • • • ^ 

Ki,X Ki,2 ^n,3 • • • ^n,n 

and B(n) = 0 . . 0 1 € J?". In the above notation, the integer n indicates the 
dimension of the state space description of the system. It is noticed that the companion 
form of a matrix defined in (2.26) is a special case of the Hessenberg form. Now consider 
the following definition : 

Definition 2.6 : A set of systems (A t (n),S(n)} will be called Hessenberg family if the 
matrices {A,} are m Hessenberg form. 
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Next we show that, if each matrix A l (n),i — 1 , 2 , ... ,1 is in Hessenberg form along with 
the same matrix B{n ) in (2.27) then they are SQ stabilizable. This class of system is 
also shown to be larger than the controllable uncertain systems which satisfy the matching 
condition. 

Theorem 2.5 A finite set of n dimensional systems {Afn). B(n)}, i = 1, 2, . . . , l m Hes- 
senberg form is simultaneously quadratically stabilizable. 

Proof : We prove the theorem by the method of induction. The matrices used in the 
following steps are assumed to be compatible in dimension. 

Step-1 : Consider a set of two dimensional system (A t (2), B(2)} in Hessenberg form de- 
fined as follows : 

1 r 

0 

1 


A, (2) = 


h \ i 1 

^■21 ^22 


and B( 2) 


Using theorem 2.1, the set of systems {A,( 2), J3(2)} will be SQ stabilizable iff for any scalar 
F and a Wi > 0 the following inequalities hold 

' (hn + F) ^ + W 1 (fi* n + F) * ] < 0 V L (2.29) 

To satisfy the above inequalities, there always exists an F such that 

(fi l u + F < 0) V i. 

By choosing such an F, one can calculate W\ > 0 from (2.29). This proves the existence 
conditions for the SQ stabilization of (A t (2), B{2)}. The controller gain is constructed 
based on the steps followed in theorem 2.1 and as given below. 


Choose Wf = FWi and W 3 such that W 3 — 1 W 2 > 0 which ensures that 

F W l W 2 1 

W = is positive definite matrix. Such a choice of U 3 is clearly always possible. 

W 2 * W 3 

Using this W matrix, there exists q such that the following holds for all the systems (using 

( 2 . 11 )) : 

A,(2)W + \V(A,(2))' - 7 B(2)(B(2))* < 0 Vi. (2.30) 
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So the required controller gain will be K = -(7/2)(5(2))*PV -1 . Here we denote this 
controller gain and the Lyapunov matrix PH, to be used in the next step, as follows : 

K{ 2) = K and W( 2) = PH. 

From the above discussions, it is concluded that any finite collection of two dimensional 
systems (A t (2),5(2)} in Hessenberg form is SQ stabilizable. This result will be subse- 
quently used in the following step. 


Step-2 : Consider a set of three dimensional system (A t (3), B( 3)} defined by : 


M 3) 



0 

1 


K 


31 


M h l 
U Z2 u 33 


A(2) 5(2) 
^21 ^22 


and 5(3) = 


0 

0 

1 


In the above partition, matrices are A z ( 2) = 


h 

hi 


li 


, 5(2) = 


0 

1 


A 1 
^21 


G 5 lx2 and A 22 G R Using theorem 2 1, the set of systems (A t (3),5(3)} will be SQ 
stabilizable iff for any matrix F and W x > 0 the following inequalities hold 


[ (A, (2) + 5(2)5) Wi + Wi (A t (2) + 5(2)5)* ]<0Vt. (2.31) 

It is noticed that the systems {Aj(2), 5(2)} are in Hessenberg form. So from the results of 
step-1, there exists an 5 and PH L > 0 such that (2.31) holds. At this point, the matrices 
5 and W\ are taken from step-1 and are given by : 


5 = K{2) and W x = W{ 2) 

Now choose PH 2 * = 5PPi and PU 3 such that PP 3 - W 2 W{ l W 2 > 0 which ensures that 


PH = 


is positive definite matrix. Using this PH matrix, there exists 7 such 


W r PU 2 
__ PU 2 * PH 3 _ 

that the following holds for all the systems using (2.11) : 


A 2 (3)PH + PH(A t (3)*) - 7 5(3)(5(3))* <OVi. (2.32) 

So the required controller gain will be K = -( 7 /2)(5(3))*PP r - 1 From the above discussion, 
it is concluded that any finite collection of three dimensional systems (A,(3).5(3)} in 
Hessenberg form are SQ stabilizable. 



2.7 Stabilization of Hessenberg Family 


30 


Step-3 : It is now assumed that a set of (n — 1) dimensional systems 

{A t (n — l),B(n — 1)} in Hessenberg family is SQ stabilizable. From theorem 2.1, it follows 

therefore that there exists an F and W x = W* > 0 such that 

[ (A,(n - 1) + B{n - 1 )F) W x + W x (A*(n - 1) + B(n - 1 )F)* ] < 0 V i. (2.33) 


The problem in hand is to show that a set of n dimensional systems {Aj(n), I3(n)} is also 
SQ stabilizable. Towards that end, we partition the special structure of A x (n) matrices of 
(2.28) as 


At{n) = 


where B(n — 1) 


h l 

r hi 

l 

0 

.. 0 




h i 
a 21 

h z 

n 22 

1 

.. 0 


A t (n - 

1) B{n- 1) " 

h l 

a n- 1,1 

Ki- 1,2 

h l 

n n~ 1,3 • 

.. i 


A 1 

L A 21 

A 1 

a 22 J 

h l 

a n, 1 

Ki 

h 1 

n n,3 

•• K,n_ 




= [° 

■ ii* 

G R n ~ lxl 

A l 21 € R lxn ~ l 

and A l 22 

€ R. 


With the help of theorem 2.1, it is concluded that the set of systems (Aj(n), B(n)} 
is SQ stabilizable if (2.33) holds. The required controller gam for the SQ stabilization of 
(A t (n), B{n)} can be arrived at using similar arguments as stated in step-1 and step-2. 
This completes the proof of the theorem. □ 


Next we shall state the stabilizability results for a set of systems {A l {n), B{n)} which 
are simultaneously transformable into Hessenberg form. 

Corollary 2.2 Consider a set of systems {Aj(n), 13(n)} having the property that there 
exists a nonsingular matrix T for which {T~ l A l {n)T,T~ 1 B(n)} are m Hessenberg form, 
then the family {A z {n),B(n)} is simultaneously quadratically stabilizable. 

Proof : The proof is similar to corollary 2.1. □ 

Remarks 2.5 : It is noticed from (2.26) and (2.28) that the controllable companion 
form is a special case of Hessenberg form. So the class of systems which are simultaneously 
transformable into Hessenberg form is larger than the controllable matched uncertain sys- 
tems. As an example consider a set of two systems {A,, B} which are not m companion 
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form defined as follows : 


x = A t x + Bu, % = 1, 2 where Ai 


0 1 
0 0 




Notice that the perturbation matrix (Ao — A L ) does not belong to the range space of 
matrix B and hence does not satisfy matching condition. But these two systems are in 
Hessenberg form. So they are SQ stabilizable by theorem 2.5. 

Remarks 2.6 : The proof of the theorem 2.5 also gives a construction procedure to 
compute the controller gain as well as the Lyapunov matrix for the SQ stabilization of 
finite number of systems {A n B} which belong to Hessenberg family. 

Remarks 2.7 : The results of this section can also be extended to multi- input systems 
using block Hessenberg form. A set of systems {A^B} in block Hessenberg form (see 
[67, 41]) is described as follows : 

h l u h\2 0 .0 

h l 2i h\ 2 h 2 , 3 • 0 

Ax = : : : : 

fyfe-1,1 ^fc-2,2 ^fc-i,3 ■ ■ h>k-i,k 

h l kl h\ 2 Ki-2,n-2 ■ ■ frfc.fc 

block matrices h hJ+ 1 are of dimension t 3 x t J+i , j 6 {1, 2, . . . , k - 1} and have full row rank 
t r In this case, t\ <tz < . . < tk = m and t x + t 2 + ■ ■ ■ + t k = n. Similar to the single 
input cases, it can be proved by the method of induction that a set of systems {A^B} in 
block Hessenberg form are SQ stabilizable. These details are omitted. 



2.8 Results For General Input Matrix Case 

In the previous section, we discussed the SQ stabilizability conditions of {A„B} by as- 
suming that the input matrix B is fixed for all the systems. In this section, we shall 
partially remove this restriction. Consider a set of systems {A l ,B t },z = 1, 2, . . . J which 
are described by the state equation / 

x = A t x + B x u (2-34) 

In this thesis, the following assumptions are made for the input matrices of the systems 

{A t ,R}. 
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Assumption 2.3 : Each of the input matrices {£,} are of the form B t = BG l where 
nonsingular matrix G t is such that (G, + G*) > 0 Vi. 

The matrix B is also assumed as 



Assumption 2.4 : The structure of static state feedback controller gain for SQ stabiliza- 
tion problem is assumed as 

K = (7/2 )B*W~ 1 (2.35) 

where 7 is a positive scalar and W is a positive definite matrix. 

The justification of these assumptions follows on similar lines as that for assumptions 
(2. 1-2. 2) (page 16) using Finsler’s lemma and the SQ stabilization of {A> B}. It is noticed 
that the input matrices {5Q are restricted 6 here in the sense that they represent particular 
types of perturbations. However this class of matrices arises in certain practical applications 
which are considered in chapter 6. 

Based on the above assumptions 2. 3-2. 4, the following result is stated for the stabiliza- 
tion of a set of systems {A l ,B l } This result is similar to lemma 2.3 and the simplified 
version of the definition 2.4. 

Lemma 2.4 A set of systems (A) Bf) are simultaneously quadratically stabihzable iff there 
exists a scalar 7 > 0 and a positive definite matrix W > 0 such that 

A t W + W A* - (7/2 )B(G t + G*)B * < 0 V 1 (2.36) 

Proof : By definition, a set of systems (A, B t } is SQ stabilizable if there exists a controller 
gain matrix K and a W > 0 such that 

(A - BJ<)*W- 1 + W- 1 (A - B t K) < 0. 

Now using K — ff/2)B*W~ 1 and multiplying both sides by the symmetric matrix W, we 
get 

A*W + WA* t - (7/2 ){BB; + £,£*) < 0 
^ A*W + WA* — ( 7 /2 )jB(G* -I- Gf)B* < 0 □ 

6 This type of uncert ain ty in input matrix also arises for matching condition with time varying uncer- 
tainty representation [59] and is discussed in chapter 1 (page 5). However, in the present case of finite 
family of systems, we are not imposing any norm type conditions. 
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Next we further simplify the above lemma in the form of following theorem. 

Theorem 2.6 A finite set of systems {A u Bf\fi — 1,2, ...,l is simultaneously quadrat- 
ically stabilizable iff there exists a matrix W x — W* > 0 and a matrix F such that the 
following inequalities hold 

(All + A\ 2 F)W x + WffAl, + A\ 2 F)* < 0 Vi 

Proof : The proof is similar to the theorem 2.1. The details are however given in Ap- 
pendix B. □ 

Remarks 2.8 : Note that in theorems 2.1 and 2.6, the existence conditions for the 
SQ stabilization of a set of systems {A t ,B t } are same as that for {A t ,B} with fixed B. 
However the controller gain will be different for the two cases. 

In the light of the above remark, the SQ stabilization of (A*, B x ] can be stated for matched 
systems similar to theorems 2.3 and 2 5. This is discussed below in brief. 

Definition 2.7 : A collection of LTI systems {A,, B t }, with assumptions 2.3, will be called 
matched systems if the submatrices A l n and A\ 2 are same for all the systems. 

Theorem 2.7 A finite set of matched systems {A l ,B l },i = 1, 2, ...,l is simultaneously 
quadratically stabilizable 

Proof : Using theorem 2.6, the proof of this result is similar to that of theorem 2.3 □ 
Similar results can be stated for the Hessenberg family of systems also. 

2.9 Controller Design Algorithm 

Based on theorem 2 1 and theorem 2.6, we present successive steps to compute static state 
feedback controller gain for the finite set of systems {A*, B t }, i = 1 , 2, . . . , l defined in (2.5 ). 
In this case it is assumed that the family of systems is SQ stabilizable. So from theorems 
2.5 -2.6, there exists a positive definite matrix W x and F such that 

(An + A\ 2 F ) W x + W x (a* u + A' 12 F)* < 0 Vj. 
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Step-0 : This step is for the computation of F and W x such that above inequalities hold. 
To compute these matrices, the following should be done 

• Compute a matrix F such that for each i, the matrix (A l u -I- A\ 2 F) is Hurwitz. 

• Compute W x by solving the Lyapunov equations 

On + A* 12 F) W x + W 1 (A\ x 4- A\ 2 F)* < 0 V *. 

In this chapter, we discuss the existence conditions of SQ stabilization for matched and 
Hessenberg family of systems. Hence, using the proof of theorems 2.3, 2.6 and 2.5, the 
matrices F and W x are obtained for these class of systems. In the next chapter, we shall 
discuss the computation of matrices F and W x for the systems referred as (a) partially 
commutative and normal systems and (b) more general family of systems. The matrices 
F and W x are used in the next step of algorithm. 

Step-1 : Using the matrices W x and F, compute Q l 0 such that 

Q l 0 = (A’ n + A\ 2 F) W x + W\ (.4* u + A\ 2 F)* for each i. Here Q l 0 is (n - m x n - m) dimen- 
sional negative definite matrix. 


Step-2 : Compute W 2 = FW X 


Step-3 : Select W 3 such that W 3 > W 2 W x 1 W 2 to ensure W = 


be a positive 


Wx W 2 
W 2 * W3 

definite matrix. Choose W 3 —r + W 2 W x l W 2 where r is a positive definite matrix. 


Step-4 : Using (2.36), compute the value of 7, and the required control law is given by 
u = -Kx = -{i/ 2 )B*W- 1 . 

The computation of 7, to satisfy (2.36), is simplified as follows. 


Define <£>, = A t W + WA; - (7/2 )B(G t + G^B * 


Q* Q 2 i — (l y/2)/m(G, + G*) 


where 


cji = A\ x W 2 + A\ 2 W 3 + HU ( A21 j + W 2(4.22 ) and 
q 2l = AS, x W 2 + A l 22 W 3 + W 2 *(A' 21 )* + W 3 {A l T2 )*. 
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Now the problem is to find a scalar 7 such that the matrices {<!>,} will be negative 
definite for each 1 — 1,2, Note that the matrix Q l 0 is a negative definite matrix 

defined in step -1. So matrices {$,} will be negative definite if [38] 

92. - (7/2 )(Gi + GO - dKOir'a < 0 V i. 

For each index i, compute a scalar 7, such that the following condition holds : 

7. > 2 A u [(92. - tiKir's-m + g.*)- 1 ] 

The required 7 is selected as the maximum of all the scalars {7,}, i.e., 

7 = max{7,}. 


Example 2.2 : In this example, a single stabilizing feedback controller is designed for a 
set of three systems {A t ,B},i = 1.2,3 using the algorithm given in subsection 2.2. The 
system parameters are given below 



' 0 

1 0 ' 


0 

1 

0 


0 

1 

0 

A\ — 

1 

-1 0 

? M — 

2.375 

-0.25 

-0.875 

1 ^3 — 

1.875 

-1.25 

0.625 


_ 0 

-1 2 


2.75 

0.5 

0.25 

1 

1.75 

-1.5 

3.25 


-2 10 

B = [0 1 2 ]*. The matrix T for the system (A 0 , B ) is given by T — 0 —2 1 

-2 12 

Using matrix T, one can transform all three systems into CC form. In the algorithm we 

’ -1 0 1 

use / = [ -2 - 3], Qo = and r = 40. 

L [ 0-1 

The matrices W and I < , computed using 7 = 3000, are given below 
1 -0.5 -0.5 

1 

and I\ = [ 75 112.5 37.5 ] . The controller gain to stabilize 


W = I -0.5 0.5 -0.5 

-0 5 -0.5 42.5 

the above three systems is I\ 0 — AT -1 


-93.75 - 75 56 25 7 


2.10 Conclusion 

In this chapter, the simultaneous quadratic stabilization problem is formulated as the 
problem of solving a set of matrix inequalities (MI). If these MI are successively solved 
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then their solution can be used to obtain a state feedback stabilizing controller for a set of 
LTI systems. Towards this end, a set of system {A z , B} with fixed B is considered. Then an 
attempt is made to identify different class of systems for which SQ stabilization problem 
is solvable, i.e., the required matrix inequalities have a solution. Following results are 
derived in this direction. The existence condition of matched uncertain system is discussed 
in the context of SQ stabilization problem. Here the matching condition is reinterpreted 
in terms of the submatrices of system matrices {A,}. This has been done by exploiting the 
structure of the input matrix B. As a special case of the matched uncertain system, it is 
shown that the class of system is SQ stabilizable if {A,} are simultaneously transformable 
into Controllable Companion form. We also characterize a class of systems which are 
simultaneously transformable into companion form. In continuation of the above results, 
it has also been proved that a collection of LTI systems is stabilizable by a single controller 
if they are simultaneously transformable into Hessenberg form. This class of systems 
is shown to be larger than the class of uncertain systems which satisfies the matching 
condition. So the simultaneous stabilizability results of this chapter (derived for finite 
number of LTI systems) are different from the existing robust stabilizability property of 
the matched uncertain systems. To accommodate the uncertainty of input matrix, we 
consider the SQ stabilization of {A,, £,}. Then under some technical assumptions on B t , 
the existence conditions of the set {A U B Z } are same as those for the set {A,,B}. The 
controller design algorithm for this finite SQ stabilizable family of systems {A,, B,} is also 
given. 

The companion form or Hessenberg form of a matrix has been used in the various design 
algorithms in control theory such as pole placement, model order reduction etc [41]. In this 
chapter, the special structure of companion and Hessenberg matrices has been explored to 
partially solve the SQ stabilization problem by static state feedback controller. 



Chapter 3 

Stabilization of Partially 
Commutative and Normal Systems 


In this chapter, two new classes of systems are introduced for which the simultaneous 
quadratic stabilization problem is solvable. The existence conditions and the controller 
design algorithm are proposed for this purpose. Any finite number of systems {A,, B t },i = 

1,2 ,1 of the form defined (2.5). are SQ stabilizable if the matrices {A* u } are Hurwitz 

and commutative or normal In addition to these systems, one more general class of systems 
is introduced here which are stabilizable by a single feedback controller. These classes 
of systems are different from the matched uncertain systems. Hence, the simultaneous 
stabilizability results of this chapter are quite different from that of chapter 2. 

The chapter is outlined as follows In section 3.1 the SQ stabilization problem of 
commutative and normal family is discussed In section 3.2, the controller design for a 
more general class of systems is derived. 


3.1 SQ Stabilization of Commutative and Normal 
Systems 

Consider the following family of LTI systems {A,, B t },i = 1,2, ... ,1 which are described 
by the state space equation 


x = A,x + B,u 


(3-1) 



3.1 SQ Stabilization of Commutative and Normal Systems 


38 


where A t G R nXn and Bi G R nXm . As in the chapter 2 (page 32), it is assumed that each 

0 


matrix B 2 satisfies assumptions 2.3, i.e., B t = 
The matrices A t are partitioned as 


1 22 



4i 


^21 

pym xm 



/, 

A 


G t and (G, + G*) > 0. 


12 


1 22 


(3.2) 


Next two lemmas show the existence of common Lyapunov matrix for a set of Hurwitz 
matrices. These results will be successively used to derive the existence conditions of the 
SQ stabilization problem 


Lemma 3.1 If a set of matrices = 1,2,...,/ are Hurwitz and commutative pair 

wise then there exists a common Lyapunov matrix P such that the following holds : 


M,P + PM t * < 0 Vi = 1,2,...,/. 


Proof : A proof of this result is given in [46] We are reproducing this proof from 
[46] for the sake of comparison For the ease of presentation, first consider two Hurwitz 
commutative matrices M\ and A/> Suppose P 0 is a positive definite matrix. Let Pi and 
P -2 be the unique positive definite solutions of the Lyapunov equations 


M\ Pi + Pi Mr = -Pq (3-3) 

il/ 2 P 2 + P 2 M; = -P v (3.4) 

Substituting for Pi from (3.4) into (3.3) and using commutativity of Mi and M 2 , we get, 

Mi(M 2 P 2 + P 2 M 2 *) + (M 2 P 2 + P 2 M 2 *)M* = P 0 

=> M 2 MiP 2 + MP 2 M 2 * + M 2 P 2 M* + P 2 M*M 2 * =P 0 
=£■ M 2 (MiP 2 + P 2 M*) + (MiP 2 + P 2 M]*)M 2 = Pq 


Since Mo is a Hurwitz matrix and P 0 is a positive definite matrix, we conclude from the 
last step that 

Mi P-2 + PoM* < 0. 

So P> is a common Lyapunov matrix for Mi and M 2 . It is also noticed that to compute 
P> we need to solve two Lyapunov equations (3.3) and (3.4). 
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Next we show that if a set of matrices {M t } are Hurwitz 1 and pairwise commutative 
then there exists a common Lyapunov matrix. Suppose given a positive definite matrix 
P 0 , let P U P 2 ,.. , P; be the unique symmetric positive definite solutions to the Lyapunov 
equations 

M t P % + P t M t * = t = 1, 2, . . . , 1. (3.5) 

Then we shall prove that the matrix Pi satisfies following Lyapunov equations 

M x P t + PiM* < 0, i = 1, 2 (3.6) 

Towards this end, define P v = M t P 3 + P 0 M*. If we can show that P v <0 for i = 
1)2, . ,1, j = i,i + 1, . . , Z, then the result follows by choosing j — l for each i. Hence, 
let i 6 {1,2,...,/}. From (3 5), P n — —P t -\ < 0. Now assume that P tJ < 0 for some 
j = (1, 2, ...,/ — 1}. Then using (3.5) we find 

M J+l P h]+l + Pt,j+iM J+1 

= Mj + i(M t Pj+i + Pj+iM?) + (M t P J+1 + p j+1 m;)m; +1 
= M t (Mj+\P J+ i + p j+1 m; +1 ) + (m 3+1 p, +1 + p j+1 m; +1 )m; 

= -{M t P J + P 3 M*) 

= —Pi] 

> 0 

This implies P lJ+ i < 0 and hence proves the results by induction. 

Note that one has to solve l Lyapunov equations to find a common Lyapunov matrix 
P Next we show that such an approach is not necessary and show a direct computation 
for a common Lyapunov matrix. 

Lemma 3.2 A finite set of Hurwitz diagonalizable commutative matrices {M t } share a 
common Lyapunov matrix and is given by p-^p -1 )* where T is a nonsingular matrix 
which diagonalizes each matrix {M t }. 

Proof : It is known from the matrix theory that diagonalizable 2 commutative matrices 
are simultaneously diagonalizable [30, 31]. Let P be the transformation matrix which 

U A set of matrices {M} are Hurwitz' means that every matrix of the^iven set {Af,} is Hurwitz. 
-Suppose matrix Ai is diagonalizable by matrix T and O i = T*' ! .4i T~* Then all nondiagonal elements 
of the matrix 0\ are zero 
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diagonalizes each matrix M t , i.e., 

6, = T~ M,T~' 

where 5, is the diagonal matrix whose diagonal elements correspond to the eigenvalues of 
matrix M % . Define P = T -1 (T -1 )*. Using this matrix the following simplification is made. 

M X P + PM* 

= T~ 1 5 l TP + PT*6*(T~ 1 )* 

= T~ l [ 5JPT* + TPT*6* } (T -1 )* 

= T~ l [ <5j + 5* ] (T -1 )* 

< 0 Vi 

The last step is true by noticing that the the matrix (<5 t + <5*) < 0 since the matrix M, is 
Hurwitz. □ 

Remarks 3.1 : In the above lemma 3.2, we use the property of commutative matrices 
that they are simultaneously diagonalizable by matrix T. Here we shall discuss how to 
compute matrix T. For this purpose, we assume that a set of commutative matrices 
{Ai,A 2 , - - , Ai} is given and one of them is diagonalizable by matrix T For simplicity, 
suppose Ai is diagonalizable by matrix T. This means that the matrix T is constructed by 
using the eigen vectors of matrix A\ and Si = T~ l A{T. Next we show that the matrix T 
diagonalizes all the matrices For this purpose, define O z = T l A t T for i = 1, 2, . . . , l. It is 
noted that the matrices {0,} are also commutative and the matrix O x = 5 X is a diagonal 

matrix. 

OiO x = 0,01 for * = 2 , 3, . . . I 
=>6 l O t = O t S l for i = 2,3,.. .1. 

The above equation ensures that the matrices O , are diagonal matrices. So the common 

Lyapunov matrix is given by P — T -1 (T -1 )*. 

Next, we show the existence of a common Lyapunov matrix for a set of Hurwitz normal 

matrices. 

Lemma 3.3 If a set of matrices {. M ,} are Hurwitz and normal then there exists a common 
Lyapunov matrix P such that the following holds : 


M t P + PM*< OVt. 
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Proof : By assumptions, the matrices {Af,} are normal, i.e., M t M* = M*M t [30], For 
these matrices, there exists orthonormal matrices T t (T~ l = T*) such that 

T t t4,r*= A,. 

Here the matrix A. t is the diagonal matrix whose diagonal elements contain the eigenvalues 
of matrix M, since each matrix M t is Hurwitz, (A, + A*) <0. 

Using P = / n , the following is obtained 

M t P + PM * 

= T* A t T, + T*A*T t 
= T*(A, + A*)T t 

< 0 

So the identity matrix is a common Lyapunov matrix for normal family. This completes 
the proof of the lemma. □ 

Next we introduce the following definition for partially commutative and normal control 
systems. 

Definition 3.1 : A set of systems {A t ,J3j will be called partially commutative if the 
submatrices {A l u }, defined in (3.2), are Hurwitz and pairwise commutative, i.e., 

^li^ii = -^n-^n ^ hi- 

Definition 3.2 : A set of systems {.4,, B,} will be called partially normal if the submatrices 
{.4^}, defined (3.2), are Hurwitz and normal, i.e., 

4l(.4n)' = (^ur^'n V *• 

Based on the above definition, we shall state the following results for SQ stabilization. 

Theorem 3.1 A finite set of systems {A„ £.}, i = 1, 2, . . . , l of the form (3.1) are simul- 
taneously quadratically stabihzable if the matrices {A^} are Hurwitz and commutative. 

Proof : From theorem 2.6, the set of systems {A t , are SQ stabilizable iff there exists 
matrices F and Hfi such that the following holds 

[ {A'n + A \ 2 F) Wi + W x (AU + A\ 2 F)* ] < 0 V * 
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By hypotheses, the matrices {A z u } are Hurwitz and pairwise commutative. So choosing 
F = 0, the above inequalities are reduced to 


Al^ + W^A^Y < OVi 

But using lemma 3.1 or lemma 3.2, there exists a single matrix Wi, such that the above 
inequality holds for all i. The construction of such a Wi matrix is discussed in the proof of 
lemmas 3. 1-3.2. Using these matrices F (= 0) and Wi in the design algorithm (described 
in section 2.9, page 33), one can compute state feedback controller □ 


Example 3.1 : Consider a set of three systems {A,, B} described in (3.1) where 


A i 


Ao — 


A, 


B 


-1 1 0 
0.5 -2 1 
1 0 1 
1 

-1.5 
0 
1 

-2.5 
1 

0 " 

0 
1 


-0.5 

0.5 

1 

-1.5 

0.5 

0 


r ' 

0.51 

1 

1.3 

0 

- 0.1 

0.5 


4 1 



4 1 

n 21 

^■22 




^12 


. ^21 

^22 . 


'Ah 

" 4-12 


. ^21 

A -22 


It is noticed from the above that the submatrices {A l u } are not same for all the sys- 
tems. So these three systems do not belong to the matched family of systems. However 
the matrices = 1,2,3 are Hurwitz and commutative. Hence these systems are 

stabilizable by a single controller. 

Using controller algorithm (described in section 2.9), we have the Lyapunov matrix 
' 0.9553 0.2277 0.0 1 

IV = 0.2277 0.7053 0.0 , and controller gain is given by K =[ 0.0 0.0 25 ]. 

0.0 0.0 2.0 _ 

Next we show the existence of a stabilizing controller for a family of partially normal 


systems. 
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Theorem 3.2 A finite set of partially normal systems {A l , B t }, i = 1,2, ... ,1 of the form 
(3.1) is simultaneously quadratically stabihzable. 

Proof : From theorem 2.6, the set of system {A l ,B l } is SQ stabilizable iff there exists 
matrices F and W x such that the following holds 

[ (A' u + A* 12 f) W, + W, (A|, + A\,F)' } < 0 V t 

By hypotheses, the matrices {A l u } are Hurwitz and normal. So choosing F — 0, the above 
inequalities are reduced to 

[A^Wx + W^Altf] <0Vj 

But using lemma 3.3, there exists a single matrix W x — / n _ m , such that the above inequal- 
ities hold for all i = 1 , 2 , . . . , l. Using these matrices F and W x in the design algorithm 
(described section 2.9), one can compute state feedback controller. This completes the 
proof. □ 


Example 3.2 : Consider a set of three systems {Ai,B} described in (3.1) where 



’ -2 

-1 

1 " 


' -1 

0 -1 " 

Ai — 

1 

-2 

0 

, a 2 = 

0 

-2 1 


1 

0 

0 


_ 0.5 

-1 1 

r- 

1 1 

0 

1 

r°i 



A3 = 


,£? = 


0 

1 


VF = 
that t 


1 -2 1 

1 1 0.5 

Using controller a: 

1 0.0 0.0 " 

0.0 1 0.0 
0.0 0.0 2.0 J 

tiese three systems are different from matched uncertain systems. 


gorithm (described in section 2.9), we have the Lyapunov matrix 
, and controller gain is given by K = [ 0.0 0.0 5 ]- It is noticed 


3.2 SQ Stabilization for More General Class of Sys- 
terns 

In this section, we shall derive results for some general class of systems. For the set of 
systems {A t , B t } the following assumptions are made 
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Assumption 3.1 : The number of inputs m > n / 2 if n is even and 
m > (n + l)/2 if n is odd. 

Assumption 3.2 : Each submatrix A\ 2 € R n ~ mxm iS a full rank matrix, i.e., 
rank (A\ 2 ) = n — m. 

Definition 3.3 : A set of full rank matrices {A\ 2 } of the partition (3.2) will be called as 
quasi sign invariant matrices if for any matrix E € A mXn_m the following conditions hold 

(A\ 2 E + (A\ 2 E)*) > 0 Vi 

As an example, suppose all the matrices A\ 2 are constant and A\ 2 = An for all i (say). In 
this case, the matrix E can be taken as E = A* u . 

Theorem 3.3 A finite set of systems {A lx B t },i = 1, 2, . . . , l is simultaneously quadrati- 
cally stabilizable if the matrices {A\ 2 } are quasi sign invariant matrices. 

Proof : For the stabilizability of the given set of systems {A t , B l }, the following inequalities 
hold from theorem 2.6 

' (a- u + a\ 2 f) w, + w, (a* u + a; 2 f)‘ ] < 0 V t. (3.7) 

Since the matrices {A' 12 } are quasi sign invariant, so by definition 3.3, there exists a 
matrix E such that 

(A\ 2 E + (A\ 2 E)*) > 0 Vi 

Using \\\ = / n _ m and F = -aE where a > 0, the inequalities (3.7) can be written as 

•^11 + (-^n)* — ot(U t + Ui) < 0 \/ i 

where U x = A\ 2 E. Since the matrix (U x + Uf) > 0 , the inequalities (3.7) will be satisfied 
by taking sufficiently large scalar a. The value of such a a can be selected as follows : 

For each index 1 , define a x = \m[ (U x + f7*) _1 (A l u + (A 5 1X )*) ] and 

a = maxlo:!}. 

Using this \\\ ( = I n _ m ) and F = -aE in the main algorithm of chapter 2 (section 2.9, 
page 33 ), one can compute controller gain. □ 
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3.3 Stabilization of a Continuously Parameterized 
Family of Systems 


In the previous sections, we have considered the SQ stabilization of finite number LTI 
systems described by 

x = A t x + B x u , t = l,2,...,h 

Here each of the above systems (A u B t ) is referred as an extreme system. Now suppose a 
finite set of extreme system {A t , Bf) is SQ stabihzable. Then the question is asked whether 
a set of systems (constructed by the positive linear combination of the given extreme plants) 
is also SQ stabihzable. 

For this purpose, define a parameterized family of systems described by 

x = A(a)x + B(a)u (3-8) 


where A{a) = ZLi a t A t and B{a) = EU a t >0Vi, and £Ui a, > 0. 
Here each matrix B t satisfies the assumption 2.3 (page 32), i.e., 


B t = BGi where B = 


0 


and (G t + G*) > 0. 


Next we prove that this parameterized family is also SQ stabihzable. 


Theorem 3.4 The parametrized, family (An B a ) described in (3.8) is simultaneously quadrat- 
ically stabihzable by a single state feedback controller if any one of the following holds 

(a) are matched systems 
or 

(b) {A l ,B l } are partially commutative or partially normal 

or 

(c) {A t ,B t } are in Hessenberg family. 

Proof : The proof of this theorem is done by employing the results of chapter 2 and 3. 
Here we are providing the proof of part (a) and the other proofs will be similar, 
proof of (a) : Since the systems {A, B t } satisfy the matching condition, there exists a 
controller gain I\ and matrix W > 0 such that the following inequalities holds by theorems 

2.7 and 2.3 

(A - B t K)*W- 1 + IT -1 (A - B t K)* < 0 V i = 1,2, . . . ,/. (3.9) 
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Next we shall show that the following inequality holds true. 

(A a - B a K)*W~ l + W~ l {A a - B a K) < 0 

Now we shall simplify the above using (3.9) in the following way. 

(A* - B a K)*W~ l + W~ 1 (A a - B a K ) 

= '£a t [{A l -B l K)*W- l + W-\A l -B 1 K ) *] 

< 0 . 

The third step is true due to (3.9) and noting the fact that positively weighted sum of 
negative definite matrices are negative definite. This ensures that the closed loop matrix 
(A a - B a K) is Hurwitz using lemma 2.1. □. 


3.4 Conclusion 

In this chapter, we solve the SQ stabilization problem of partial commutative family and 
partial normal family of LTI systems. The existence conditions are ensured for these 
class of systems. The results for some more general class of systems are also derived. It is 
important to note that these classes of systems are different from that of matched uncertain 
systems. Hypothetical numerical examples are also provided to verify these results. 



Chapter 4 

Simultaneous Stabilization by 
Output Feedback 


In this chapter, new classes of LTI systems are identified for the s im ult an eous quadratic 
stabilization problem by static output feedback. 


4.1 Introduction 

Design of static state feedback controller for the control system design has been discussed in 
the previous two chapters. In these cases the feedback controller of the form u = — Kx(t ) 
has been used. This type of controller requires knowledge of only instantaneous value 
of state x(t) and hence one must have all the state variables x(t) available for feedback. 
However, in practice, it is not always possible to measure all the states. So the design of 
static output feedback controller by using measurable states is a more realistic approach 
in control theory. 

Recently a sufficient condition for output feedback stabilizability of linear matched 
uncertain systems has been developed in [65, 64, 83, 27]. It has been established that 
if the nominal system is minimum phase and has relative degree of unity, then robust 
stabilization is possible using static output feedback controller [65, 83]. In [27], this result 
is extended to multivariable systems with equal number of inputs and outputs. In this 
chapter, we report new classes of systems which will be SQ stabilizable by static output 
feedback. 

Consider a family of multi-input multi-output (MIMO) systems {Ai,B u C} * = 1 , 2 ,...,! 
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described by the state equations, 


x = A t x + B z u, y-Cx (4.1) 

where x € R n is the state vector, u e R m denotes control input and y e R p is measured 
output. To stabilize the above set of systems, the following output feedback controller is 
considered 

u = —Ly[t). (4.2) 

The basic problem is to design a matrix L such that the closed loop matrix for each system 
(A — B X LC) is Hurwitz, i.e., all the eigenvalues of {A x — B Z LC) are in the open left half of 
the complex plane. 

Before going to the actual discussion, here we show how dynamic output feedback can 
conceptually be replaced by static output feedback [7]. More specifically, suppose we want 
to design a v th order dynamic output feedback controller of the form 

x c = .4 c :r c + B c y, u = C c x c + D c y. (4.3) 


The augmented system (using 4.1 and 4.3) can be written as 


X a — A al x a T Bmlla, y a — C a x a , u a Kya 


where x a = [x* £*]*, u a = [u* x* c }*, y a = [y* £*]*, 


Ani — 


A t 0 

0 0 


Bai = 


B t 0 
0 I v 


, C a 


C 0 
0 I v 


,K 


So the design of dynamic output feedback controller is equiva! 
output feedback controller for the augmented system. 


(4.4) 


D c C c 

B c A c 
ent to the design of static 


4.2 Problem Definition and Results 

In general the design of output feedback controller is difficult both from the existence con- 
dition as well as algorithmic point of view. So for the design of matrix L e R mxm , we have 
combined the Lvapunov stability approach and some mathematical tools for the solution 
of matrix inequalities. 
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Definition 4.1 : A set of systems (4.1) is simultaneously quadratically stabilizable by 
static output feedback if there exists a matrix W and matrix L such that the following 
inequalities holds 


(At B l LC)*W 1 + W — B l LC ) < 0 V i. □ 

For the set of systems (4.1), following assumptions are made. 


Assumption 4.1 : The matrix B t = BG t where nonsingular matrix G l E R mxm is such 
that {Gi + G *) > 0. 


outputs, i.e., p = m. The matrix B 


>mxm 


IS 


Assumption 4.2 : The number of control inputs is equal to the number of measurable 

0 

In 

L. 

a nonsingular matrix. These particular structures of B t and C matrices ensure that the 
determinant of the matrix [ CB r ] is not equal to zero, i.e., det[CBi] ^ 0. 


and C = [Ci C 2 ] where C 2 £ R 


Assumption 4.3 : The transfer function matrix C(sI n -A l )~ l B l for each system is strictly 
minimum phase i.e., all the zeros of the system will lie in the open left half of complex 
plane. 

In the above assumption 4.1, the particular structure of B matrix is always possible 
using a similarity transformation. Notice that these assumptions (4.2 -4.3) are made for 
the sake of convenience. The expression det[CB l ] is known as the high frequency gain of 
the system (A,,£„C) [27]. 

To design output feedback controller, we consider the state feedback controller described 
in chapter 2. In chapter 2, the state feedback controller, of the form 
u = — Kx(t ) = —~fB*W~ l x(t) where 7 is a positive scalar, is used. The W matrix is the 
positive definite solution of the following matrix inequalities (from lemma 2.4) 

A t W + WA; - ('j/2)(B l B* + BB;) < 0 V i (4.5) 

In the above (Riccati) inequalities, 7 is a free parameter which must be chosen properly 
to ensure the positive definiteness of W . When all the states are not available, one needs 
to consider the output feedback controller. Due to the formulation (4.4) above we look 
for ways of designing the output feedback from the state feedback law. In view of (4.5) 
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and the above discussions, we restate the problem of SQ stabilization by output feedback 
controller as follows : 


Definition 4.2 : A set of systems in (4.1) is simultaneously quadratically stabilizable by 
static output feedback if there exists a positive definite matrix W and a matrix L such 
that the following conditions hold for some positive scalar 7 : 


(а) A t W + WA* - (7/2) (B Z B* + BB '*) < 0 Vi (4.6) 

(б) LC = (7/2 )B*W~\ (4.7) 


In the above definition, the state feedback controller K = (7/2 )B*W~ 1 guarantees the 
stability due to (4.6), while due to (4.7) there is an output feedback gain L corresponding 
to state feedback gain K. Note that if (4. 6-4.7) hold for each i = 1,2,..., I then V ( x ) = 
x*W~ l x will be a quadratic Lyapunov function for each closed loop system. 

Let the matrices A t , B t , C be partitioned as follows : 


A. = 


A * A 1 
Ai 71 12 


A 


21 


*•22 


, B z = BG t = 


G u C = [Ci C 2 ] 


(4.8) 


where A l n € R n ~ mxn ~ m , A\ 2 <E R n ~ mxm and G z is a nonsingular matrix. 


Theorem 4.1 Suppose the system (4-1) satisfies the assumptions 4-1-4-3 and matrices 
A t , B t , C are partitioned as given in (4-3)- A set of systems {A 2 , B % , C} are simultaneously 
quadratically stabilizable by a output feedback control law u = -Ly with a nonsingular 
matrix L iff there exists a positive definite solution kVi of the following inequalities 

(4, - A\,F)Wi + If", (4, - A[ 2 FY < 0 V i (4.9) 


where F = C 2 l Ci- 
Proof : (Sufficiency) : 

By hypothesis, there exists a matrix W x = Wf > 0 such that the following inequalities 
hold, 

P l n - A l n F)W y + WfiAl, - Ai 2 F)*] < 0 V *. 

We shall prove that there exists a matrix W and a nonsingular matrix L such that in- 
equalities (4. 6-4. 7) hold. For this purpose, we construct a positive definite matrix W m 
the following way. With the help of matrices F and Wi, we choose 

w; = -FW\. 


(4.10) 



4 2 Problem Definition and Results 


51 


Then, we select W 3 such that W 3 — W 2 W X 1 W 2 > 0 which ensures that W 


Wi W 2 
\Y* W z 

is a positive definite matrix. Such a choice of W 3 is always possible as shown in the proof 
of theorem 2.1 of chapter 2. 

Using this matrix TV , consider the following matrices 


^ = A,w + vfa : - (7/2 )( B t B * + bb ;). (4.11) 

Now the problem in hand is to show that {$*}, i = 1,2, . .,1, are negative definite and 
there exists a nonsingular matrix L such that (4.7) is satisfied. 

Now from (4.10) the matrix F = — W 2 Wf 1 while by hypothesis of the theorem matrix 
F = C 2 l C\. Hence (4.10) implies : 


We partition <T> t 


'll rU2 

and using F = -W 2 \\\ l = C 2 l C u we get 


w; = -c 2 l c x w x =» c x w x + c 2 w 2 * = 0. (4.12) 

with 4 >\ j € FL n ~ mxn ~ m . Then, after simplifying (4.11) 


#1 ^ 


fti = (^U - A\ 2 F)W x + WMh - A\ 2 F)* (4.13) 

<t > l 22 = (A^w 2 + A^W 3 ) + (A' 21 W 2 + A^Wz)* - (7/2 )Im(G % + Gt) (4.14) 

d>\ 2 = A\ l W 2 + A l 12 W 3 + Wi«)* + W 2 (Ay* (4.15) 


The matrix T, will be negative definite if [38] 


(/) fa < 0 and (II) fa - (fa)* (fa)~' fa < 0. 


Observe that the matrix fa is negative definite for each i = 1 , 2 , . . . , l by the hypothesis 
of theorem. Also observe that fa and fa are independent of 7 and the matrix 
(Gi + G*) > 0. So the matrix fa can be made negative definite by choosing a sufficiently 
large value of 7. Since o* n and fa are independent of 7, we can choose 7 such that 

fa - (fa)* (4) -1 fa < 0 Vi. 

wfc'NTRAL LIBRARY 

I I. T.. < AN PUR 

••Ufa. A 123585 
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The value of 7 is computed as follows : 

Using the expression of </> 22 from (4.14), the above inequalities are rewritten as 


(7/2) [G, + G: ] > [S t - (ft 2 )* O^)- 1 <f>\ 2 ] V 2 (4.16) 


where matrix S , = (A* 21 W 2 + T? 22 IT 3 ) + (A l 2l W 2 + A\ 2 W 3 )*. 
Now for each index 2, choose a scalar 7, such that 


h > 2A m . (5 t — (^* 12 )* (<^ u ) 1 (j>\ 2 )[Gj + G*] 1 j . 


Let the value of 7 = max{7,}, 2 = 1,2 ,...,/. For this value of 7 and W, the matrix 
defined in (4.11), will be negative definite for all i. 

The existence of nonsingular matrix L is discussed as follows. 

The matrix W is a positive definite matrix. So CWC* > 0. By partitioning C and W we 
write 


>1 Qt] 


Wi 

W 2 * 


W 2 

W 3 


c* 

Co* 


> 0. 


Using (4.12), the above inequality is further simplified and written as 


C,W 2 + C 2 W 3 


c 2 * > 0. 


Since the matrix C 2 is a nonsingular, the above inequality ensures that the matrix 
[ CiW ‘2 + C 2 W 3 ] is a nonsingular matrix. 

Let us define the matrix 

L = (7/2) [ CiW 2 + G2W3 ] _1 - ( 4 - 17 ) 

Using this matrix L, we observe that (using (4.12)) 

LCW = (7/2) [ C\Wi + C2W3 ] -1 [C\Wi + G 2 1F 2 C1W2 + C2W3] 

= (7/2 )B\ 

Thus using the matrix L in (4.17), equation (4.7) is satisfied. Hence the controller 
u = - Ly(t ) stabilizes the family of systems. 


Necessary : Suppose there exists a W and a nonsingular matrix L such that (4.6) and 
(4.7) hold. Consider the matrix 

4>, = A,W + WA‘ - (7/2 + BB’) < 0. 


(4.18) 
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Now from (4.7) (note that B = 


0 


and L is a nonsingular matrix), we write 


CW = ( 7/2 )L~ 1 B* => C X W X + C 2 W 2 * = 0 , CxW 2 + C 2 W 3 = ( 7 / 2 )IT 1 (4.19) 


Using (4.19), we partitioned the matrix = 
0* , is given by 


01 1 012 
( 012 )* 022 


of (4.18) where the submatrix 


0 ii = (^11 " ^12 W + Wi(A* u - a\ 2 ft 

Since the matrices {<l> t } are negative definite, so are the matrices {q l n }. Hence the proof. □ 


Remarks 4.1 : We discussed the static output feedback stabilizability of square systems 
(equal number of inputs and outputs). These results can be applied to nonsquare systems 
by squaring down the systems. The results of squaring down of nonsquare LTI system 
reported in [58] may be useful for this purpose. 

To verify the existence conditions given in (4.9) of theorem 4.1, one needs to solve the 
following Lyapunov equations for the existence of a single matrix W\ : 


( 4 , - A' n F) w, + w, (4, - 44)' < 0 V !. 


However, the necessary condition for the existence of a common solution of the above 
inequalities is that the matrices {A l u - A\ 2 F} must be Hurwitz for each i. Due to the 
above requirement, we shall state a lemma which will be useful to prove the existence 
conditions given in (4.9). In the following development, the notion of system zeroes [56] 
are used to define minimum phase systems. A system will be called minimum phase system 
if the zeroes of the given system lies in the open left half of complex plane. 

Lemma 4.1 Suppose that the matrices A 1 ,B l ,C are partitioned as in (4.8). Then the 

transfer function C{sl n - A t )~ 1 Bi is minimum phase iff the matrix 

(^ u _ A\ 2 C 2 1 C x ) is Hurwitz, i.e., all the eigen-values of matrix (A\ x - A\ 2 C 2 l C x ) have 

negative real parts. 


Proof : This result is based on [26] and proved as follows. 
From the definition of zeros [56], rank 


H 

I 

6 

1 

CO 

-^12 

“^21 

slm ^22 

-Ci 

— (-2 


0 

Gi 

0 


will decrease for 
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any value of s which is a zero of the system realization {A u B u C}. It is also known that 
multiplication of a nonsingular matrix does not change the rank of a matrix. 


Define the matrix 



■Sl n -m -4] j 

_ A * 

^12 

0 


m 

0 

0 

Z- 

~^21 

qT _ A} 
* 1 m SA‘22 

G, 


—C 2 l Ci 

•Ln 

0 


-Cl 

-C 3 

0 


0 

0 

cr 1 


sln-rn ~ -4* u + -4 a ^2 *C\ 
* 

0 


* 

I, 

'2 


0 


m 

Co 0 


where the submatrices denoted by * can as- 


sume any form. The rank of above matrix will decrease if s is an eigenvalue of matrix 

(-4 a — A\ 2 C 2 l Ci). □ 

In theorem 4.1, we derive existence conditions for the SQ stabilization of a set of system 
{Ai, B t ,C} by a static output feedback controller. These conditions are derived in terms 
of matrix inequalities defined in (4.9). But what are the structures of {. A t ,B t ,C } such 
that SQ stabilization problem is solvable by output feedback? The rest of this chapter is 
devoted to obtain such classes of systems. 


4.3 Stabilization of Matched Systems by Output Feed- 
back 

The concept of matching condition for linear time varying uncertain system is discussed in 
chapter 1 (page 5). Using matching condition 1 , robust stabilization by output feedback is 
solved in [27]. However, in the present context, we shall call a set of LTI system {A t , B t ,C} 
as matched system if the submatrices {A^} and {A\ 2 } are same for all the systems. The 
formal definition is as follows. 

Definition 4.3 : A set of LTI systems {A u Bi,C} will be called matched system if the 

submatrices {4 l 1 u } and {.4 l 12 } are same for all the systems. 

Defining this, the SQ stabilization results of such systems is stated as follows. 

Consider an uncertain systems x = (do + Ad)x + (Bo + A B)u,y = Cx. This system is said to satisfy 
matching condition [27] if A.4 = BqO\ and A B = BqO 2 where HO 2 II < 1- 
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Theorem 4.2 A finite set of matched systems {A x , i B x , C}, i = 1,2, ... ,1 described in (4-1 ) 
is simultaneously quadratically stabilizable by static output feedback controller. 


Proof : To prove the above theorem consider the special structure of A x and B x matrices 
as partitioned in (4.8). By definition 4.1, the matrices A x and B x defined for matched 
system are as follows 


■^11 

A \2 

r—f 

A 1 

^22 


and B, 


0 

dm. 


G.. 


For the stabilizability of the given set of systems, the following expression will be negative 
definite, i.e., 

[ (A n - A U F ) W x + W l (A u - A 12 F )* ] < 0 (4.20) 


where F = C 2 l C\. 

This is a Lyapunov inequality. Since {A l ,B l ,C} are minimum phase, the matrix 
(.4 U - Ai 2 F) is Hurwitz by lemma 4.1. Hence, there exists a matrix W\ > 0 such that 
(4.20) holds. The matrix W\ is computed by solving the following Lyapunov equation 


( (An - A n F) W x + W x [A xx - A l2 F )* ] - Q 0 


where Q 0 is any negative definite matrix of dimension (n — mxn — m). The matrices W x 
and F are used in the step 1 of the algorithm (described in section 4.5, page 57) to find 
controller gain. □ 

As a special case of theorem 4.3, it is proved in next theorem that if each matrix 
A x ,i = 1 , 2 ,...,/ is in block companion form [62], then the systems {A X ,B X , C} are SQ 
stabilizable. 


Theorem 4.3 A finite set of systems {A x , B x ,C},i = l,2,...,l in block controllable com- 
panion form is simultaneously quadratically stabilizable by a static output feedback. 


Proof : To prove the above theorem we consider the special structure of A x and B x ma- 
trices of (4.1) as given below, 


0 I m 0 ... o 

0 0 J m ... 0 

0 0 o ... d m 

—OIq ~ o) ‘ ■ * — ®r 


-An A.12 
A-21 A 22 


(4.21) 
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0 


and Bi — \ \ G t where r = n/m. 

Note that for each system the matrices An and Ayi are same. Using theorem 4.2, we 
can conclude that {.A*, B t ,C} are SQ stabilizable by an output feedback controller. This 
completes the proof of the theorem. □ 


4.4 Results for Mismatched Systems 

In the previous section, we discussed the stabilization of the matched systems by output 
feedback controller. In the present section, we identify two new classes of systems which 
are stabilizable by single output feedback controller. These classes of systems are different 
from matched systems. Towards this end, we shall prove the following results for the SQ 
stabilization problem by output feedback. 

Theorem 4.4 A finite set of LTI systems {A l , B t ,C}, as partitioned m (4-8), is simulta- 
neously quadratically stabilizable by static output feedback if the matrices {A l u — A'nC^C,} 
are Hurwitz and commutative. 

Proof : From theorem 4.1, the set of system {A u B t ,C} is SQ stabilizable iff there exists 
a positive definite matrix \V X such that the following holds 

[ (.-I’ll - A^F) Wi + W, (A' u - A' U F)' } < 0 V I (4.22) 

By the hypotheses of the theorem, the matrices {A* u - A\ 2 F} are Hurwitz and pairwise 
commutative. Hence using lemma 3.1 of chapter 3, there exists a single matrix W*, such 
that (4.22) holds for all i. This matrix W x is used in the design algorithm (described in 
section 4.5, page 57), to compute output feedback controller. □ 

Theorem 4.5 A finite set of LTI systems {A l ,B l ,C}, as partitioned in (4.8), is simulta- 
neously quadratically stabilizable by static output feedback if the matrices {A 1:L — Ci} 

are Hurwitz and normal. 

Proof : From theorem 4.1, the set of system {A i ,B 1 ,C} is SQ stabilizable by a static 
output feedback controller iff there exists a matrix W x such that the following holds 

[ (4, - .-ikcr'Ci) w, + w, (4, - 4 2 c 2 ->c,)' ] < 0 v i 


(4.23) 
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By the hypotheses of the theorem, matrices {A J U — A^C^Ci) are Hurwitz and normal. 
Hence using lemma 3.3 of chapter 3, there exists a single matrix Wi = In-m> such 
the (4 23) holds for all i. Controller gain is computed using this matrix W\ in the design 
algorithm (described in section 4 5, page 57). This completes the proof □ 

Remarks 4.2 : It is observed from definition 4.1 that a set of systems {A Z ,B Z ,C} are 
matched if the submatrices {A^} and {A* 12 } are same for all the systems. However, m 
theorem 4 3 or theorem 4.4, we consider a class of systems for which the submatrices {A l u } 
and {A l 12 } are not needed to be same for all the systems. In this sense, they are different 
from the matched systems. This will be also clear from the example 4.1. 


4.5 Controller Design Algorithm 

Based on the above proof of theorem 4.1, we present successive steps to compute a static 
output feedback controller gain for the finite set of systems {A t , B z ,C},i — 1,2,..., I 
defined in (4.1). In this case it is assumed that all the systems are SQ stabilizable by an 
output feedback controller. So from theorem 4.1, there exists a positive definite matrix W 1 

such that + 

(A 1 !! - A\ 2 F) W x + Wi (A l n - A\ 2 F) < o Vi 

where F = C 2 l Ci. 

Step 0 : This step is for the computation of W, such that above inequalities hold. In 
the previous sections of this chapter, existence conditions (hence the construction of Th 
matrix) are discussed for the following cases : 

• Matched systems 

. Matrices {A* u - A[ 2 F} are Hurwitz and pairwise commutative 

• Matrices {A l u - A\ 2 F} are Hurwitz and normal. 

The Wi is subsequently used for the controller design purpose. 
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Step-1 : Using matrices W x and F, compute QJ such that 

Q o = (^u — A\ 2 F) \\\ + W\ ( A l u — A\ 2 F)* for each i. Here Q l 0 is (n — m x n — m) dimen- 
sional negative definite matrix 


Step-2 : Compute W 2 = -FW X 


Step-3 : Select IV3 such that W3 > W 2 W x 1 W 2 to ensure W 


be a positive 


Wi W 2 
W 2 * Wz _ 

definite matrix. Choose W3 = R + W 2 Wf 1 W 2 where R is a positive definite matrix 


Step-4 : Using (4.6) and (4.7), compute the value of 7, and the required controller gain 
is given by L = {'y/2){C\W 2 + C 2 Wz)~ l . The computation of 7 to satisfy (4.6-4.7) is 
simplified as follows. 


Define <i>, = TJU + W A* - (7/2 ){B l B*+BB;) 


Q 0 

qt to-(7/2)(G. + G7) 


where 


q t = A\ j W 2 + A\ 2 W 3 + + W 2 {A l 22 )* and 

ft, = AH n W 2 + ASaWi + W' 2 W + Wz(A 22 )*. 


Now the problem is to find 7 such that the matrices {<U} will be negative for each 
» = 1 , 2, . . . , l. Note that, the matrix Q l 0 is a negative definite matrix defined in Step-1. So 
the matrices {<!>,} will be negative definite if [38] 

(t) Qo < 0 (^) Q2i ~ (7/2) (^i + G*) — q*{Qo) 1 ft<0Vi. 

Now, For each index i, compute such that the following holds . 

7, > 2Ajvf[ [ ft, - ql (Qo) -1 ft ] [ + G * ] _1 ]■ 

The required 7 is selected as the maximum of all such scalars {7i}> J - e -> 

7 = max{7i}. 


This completes the design algorithm □ 
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Example 4.1 : Consider a set of three systems {A l ,B l ,C},i = 1,2,3 where 


II 

-2 -1 0 

2 -2 1 

II 

CM 

-1 0 0 " 

1 -2 1 

ll 

CO 

’-11 0 

2 -2 1 


1 12 0 


0.5 -1 1 

L 


1 1 0.5 _ 


£ = [001]* and C = [ 1 0 1 ]. 

Here matrices {.4 l u — A^C^Ci} are Hurwitz and normal. So by theorem 4.4, there exists 
a static output feedback controller. Using the algorithm given in section 4.5, we compute 


controller gain L = 5 and the Lyapunov matrix W = 


1 0 -1 
0 1 0 
-10 3 


It is also noticed that the submatrices {A^} are not same for all three systems. So 
those systems are different from matched systems. 


4.6 Conclusion 

In this chapter, three new classes of systems are identified which are simultaneously 
quadratically stabilizable. An algorithm is also provided to compute such a controller. 
These classes of systems are different from the matched uncertain systems and considerably 
improve the results reported in the literature [27, 65]. Although these conditions are 
devised for linear systems using constant output feedback, they are equally applicable to 
linear systems using dynamic output feedback. 



Chapter 5 

Stabilization of Time Delay and 
Time Varying Systems 


This chapter is divided into two parts. In the first part, sufficient conditions for memoryless 
state feedback stabilization of a set of time delay systems are presented. The results of 
chapter 2 and 3 are extended to stabilize the delay system. Lyapunov stability approach 
is employed to investigate this problem. Furthermore, these results also allow one to 
synthesize a controller using a simple algorithm. 

In the second part, we discuss the usefulness of SQ stabilization problem for time 
varying systems. More specifically, we characterize certain classes of time varying systems 
which are stabilizable by state feedback controller. This chapter is organized as follows. In 
sections 5.1- 5.3, we discuss the stability of time delay systems. The results of time varying 
systems are discussed in section 5.4. 


5.1 Preliminaries 

In practice, time delay is commonly encountered in various engineering systems such as 
chemical processes, man-machine systems, long transmission lines m hydraulic and rolling 
mill systems etc. Hence the study of stability problem of time delay systems has drawn 
the attention of researchers over the past decade [39, 28, 82, 69, 17]. 
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Consider a set of state delay system described by 

dx(t) 

= ±(t) = A t x{t) + A h x(t — h) + Bu (5.1) 

x{t) = t ?( i),i € [-h, 0], 

where x € R n is the state vector, u € R m is the control input and A t ,B,Ah are constant 
matrices. The parameter h > 0 is a scalar number representing the delay in the system and 
7j(t) € C[— /t, 0] is initial state vector-valued function which is continuous on the segment 
—h < t < 0. We represent x t — x(t + 0) where 9 € [— h, 0]. This set of systems will 
be denoted by {.4,, B,Ah}. It is also assumed that the solution of (5.1) exists. For more 
details regarding this see [39, 28, 82]. 

Sufficient conditions are derived in [10, 16, 28, 82] to stabilize a single Time Delay (TD) 
system (.4,, B, Ah) using memoryless state feedback controller of the form u = -Kx(t). 
This type of controller requires knowledge of only instantaneous value of state x(t). Using 
Riecati equation approach 1 , the stabilization of uncertain system is solved in [11, 55, 61]. 
In this approach linear state feedback controller is sought to stabilize the delay system. The 
existence condition of these results are related to the solution of a parameter dependent 
Riccati equation. If the Riccati equation has some solution, then controller gain can 
be computed by their algorithm. It is also known that matched 2 uncertain system is 
asymptotically stabilizable by state feedback controller [61, 44]. 

In this part of the chapter, different class of delay systems are identified which will be 
stabilizable by a single memoryless state feedback controller. The state feedback controller, 
to be designed, is of the form 

u = —Kx(t). (5-2) 

In order to design matrix I< € R mxm , we have combined the Lyapunov stability approach 
and some mathematical tools for the solution of the Riccati equation. 

The following assumptions are made in this chapter about the family of systems described 
in (5.1). 


^ 1 (0)» where A* t ). 

“ * 441 

A.4(t) = B 0 Qx(t), A.4 h = BoOs(t), AD = B 0 O 3 (t) where ||Qs(*)ll < 
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Assumption 5.1 : The matrix B is assumed to be of the form B = 

dm 

Assumption 5.2 : The delay matrix Ah = BD where matrix D € R mxn . 

Assumption 5.3 : It is assumed that the system pairs {A U B} are stabilizable, i.e., the 
non delay systems are stabilizable. 

In the. assumption 5.1, the particular structure of B matrix is always possible using a 
similarity transformation and can be justified on the line of justification of assumption 2.1 
(page 16). 


5.2 Stability of Time Delay Systems 

In this subsection, we shall discuss the stabilizability of a single TD system (A l: B,Ah) 
bv state feedback controller. Using the controller u = -Kx(t), the closed loop system is 

written as 

x(t) = (4 - BK)x(t) + A h x(t - h) (5-3) 

To verify the stability of the closed loop system (5.3) (after plugging the controller), we 
consider a particular type of Lyapunov functional 

V(x t ) =x*(t)Px(t)+ f x*(s)Qx(s)ds (5-4) 

J t—h 

where P = P* > 0 and Q = Q* > 0 are positive definite matrices. Notice that V(x t ) > 0 
since P > 0 and Q > 0. Next we shall introduce the following definition regarding the 
restricted version of the stability of delay systems. 


Definition 5.1 : The closed loop TD system (5.3) is quadratically stable if there exists a 
quadratic functional V(z.) of the form (5.4) and the derivative V(* t ) is negative definite 
along the trajectory of closed loop system. 


Remarks 5.1 = For quadratically stable system, the solution x(t) of closed loop system 

converges to x(t) = 0 as time t -* oo [28, 39]. 

With this brief introduction, we shall state a lemma for the stabilrzability o a sin- 
gle TD system (.4,, fl, .4„) by state feedback controller. In this context, a delay sys e. 


5 2 Stability of Time Delay Systems 


63 


will bo stabilizablo by state feedback controller if the corresponding closed loop system is 
quadratically stable. 

Lemma 5.1 A time delay system (A t ,B,An) is memoryless state feedback stabihzable via 
the Lyapunov functional (5.4) if there exists a memoryless controller u = — Kx(t ) and a 
positive definite matrix P such that the following inequality holds : 

(.4 1 - BK)*P + P(4 - BK) + PA h Q~ l A* h P + Q< 0. 

Proof : Choose the Lyapunov functional as given in (5.4) 

V(x t ) = x*(t)Px(t) + I x*(s)Qx(s)ds. 

J t—h 

Defining A t = .4, - BI\ , the derivative of V(x t ) along the solution of closed loop system is 


V(x t ) = x*(t)(A* l P + PA l + Q)x{t) 

+x *(i - h)A* h Px(t) + x*{t)PA h x{t -h)-x*(t- h)Qx(t - h) 


= [ x*[t) x*(t - h) ] 

= [ x'(t) x*(t — /i) ] [ THT* ] 


A* h P 


A:P + PA t + Q PA h 

-Q 

x{t) 

x(t — h) 


x(t) 

x(t — h) 


Here matrix T 


In ~PA h 
0 Q 




A*P + PA l + PA h Q~ l A* h P + Q 0 
o -Q~ l 

Since matrix T is nonsingular, V(x t ) .will be negative definite if H < 0. So delay system 
will be stable if 

(A t - BI\)*P + P(A t - BK) + PA h Q~ l A* h P + Q< 0 (5.5) 


□ 

Lemma 5.1 provides a delay independent stability criteria for a single system ( A h B,A h ) 
since the above inequality does not include the delay h. 

In this chapter, we use Q = L in the general expression of V(x t ) in (5.4). So using Q = I n 
in (5.5), the following definition is introduced here. 


Definition 5.2 : A set of TD systems {.4, B, A h } is SQ stabihzable if there exists a matrix 
K and a matrix P > 0 such that the following hold 

(At - BI<) m P + P(Ai - BK) + PA h A* h P + I n < 0 V *. 
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Next wo shall state the simultaneous stabilization of a set of TD system {A u B , Ah) in 
the form of following lemma 

Lemma 5.2 A set of systems {A u B, Ah) defined m (51) is simultaneously quadratically 
stabihzablc by a state feedback controller if there exists a scalar 7 and a positive definite 
matrix W such that the following inequalities holds 

A l W + W A* - 7 BB* + A h A* h + WW <0Vi. (5.6) 

Proof : Suppose there exists a matrix W and 7 such that (5.6) holds. Now pre- and post- 
multiplying (5.6) by W~\ we get 

{A t - BK)*P + P[A t - BK ) + PA h A* h P + I n < 0 Vi (5.7) 


where P — W and 


(5.8) 


K = { 1 /2)B*W~ 1 . 

Next choose a Lyapunov functional as given in (5.4) with Q = I n which yields 

V{x t ) = x*{t)Px{t) + f x*(s)x(s)ds. 

J t—h 

Thf> derivative of V' (;;;,) will be negative definite for each system (A t ,B,Ah) in the light of 
•qualities (5.7) and lemma 5.1. So each system is stabilizable by a control law u = -Kx. 


1IK 

□ 

u 


It is observed that if (5.6) holds, then simultaneous stabilizing controller is given by 
_ _(^/2 Next the following partition has been made in this section. 


. 4 , 



^12 

4 i 

^22 


,B = 


0 


,Ah = 


0 

D 


(5.9) 


where A' n € R n ~ mxn ~"\A\ 2 € R n ~ mXin ,De R 


mxn 
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5.3 Main Results 

Using above partition (5.9), the following result is derived. 

Theorem 5.1 Suppose the system (5.1) satisfies the assumptions 5. 1-5. 3 and matrices 
{/h, B, A i , } are partitioned as given in (5.9). These system are simultaneously quadratically 
stabdizablr if there exists a matrix Wi = W* > 0 and a matrix F such that the following 
in < ‘tpi all ties hold 

(A\ 1 + + W X {A ‘ u + A\ 2 F)* + VUi(I n _ m + F*F)W, <0Vi (5.10) 

Proof ■ The proof is similar to the proof of theorem 2.1 of chapter 2 and can be found 
in Appendix B. □ 


Remarks 5.2 : In theorem 5.1, we assume the matrices B and Ah are constant for all 
the systems. However, we can derive similar results for a set of TD systems {A l ,B t ,A h i} 
where 

Bi = BG U (G t + G*) > 0 andA/u = BD t 

The details of these issues are not discussed here. 

The condition derived in the above theorem is not appealing enough because it involves 
two additional matrices \\\ and F. However, it is expected that the condition should be 
derived in terms of given state space parameters A, B, and A h . For this purpose, the 
following class of uncertain TD systems are considered for the solution of SQ stabilization 
problem. These svstems are (1) matched systems, (2) partially commutative sys- 
tems and (3) partially normal systems. Each of these family of systems satisfies the 
assumptions (5. 1-5.3). These classes of systems are defined in the context of finite number 
of TD systems as follows and similar in definition given in chapter 3. 


Definition 5.2 : A set of TD system {A„B,A h } is said to be matched system if the 
submatrices end A',, (defined in (5.9)) are same for each system. 


Definition 5.3 : A set of TD system {.4„ B A,.} is said to be partially commutative if 
the sulmiatrices {.T u } are Hnrvvitz and pairwise commutate. e. 
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Definition 5.4 . A set of TD system {A I5 B, A is said to be partially normal if the 
submatriees {A l u } are Hurwitz and normal. 


Defining the above classes of systems, we shall state three theorems for the solution SQ 
stabilization problem of these class of systems. 

Theorem 5.2 A set of matched TD systems {A l: B. A^} is simultaneously quadratically 
stalnhzahlc by a single memoryless feedback controller. 


Proof : Based on the definition 5.2, we define and partitioned the matched systems 
{A„B,A /t } as follows : 


A = 


^11 

•A 12 

An 

^22 




0 


>^ 4/1 


0 

D 


It, is noticed that the non delay systems {A t ,B} are stabilizable, so (A U ,A 12 ) is also 
stabilizuble by lemma 2.3. Thus there exists a matrix F such that (An +AuF) is Hurwitz. 
Using such a matrix P, solve for Pi = P* > 0 from the following Lyapunov equation 


(,l n + A n FyPi + P X (A U + A U F ) + (I n _ m + F*F) + Q = 0\/i 


where Q is any positive definite matrix. By pre- and post-multiplying above equation by 
matrix \\\ = Pf \ the condition given in (5.10) is satisfied. Hence the proof. □. 

Theorem 5.3 A set of partially commutative TD systems {A„J3,A/ l } is simultaneously 

quadratically stabilizable.. 

Proof : To prove this theorem, take F = 0. So (by substituting F = 0 in (5.10)) the 
problem is to find a matrix = W* > 0 such that the following inequalities hold 

+ Wl(4l)‘ + WVn-mW, < 0 V i. (5.11) 

By hypothesis, the matrices {An } are Hurwitz and pairwise commutative. Hence by lemma 
3.1 or 3.2 of chapter 3 (page 39), there exists a matrix Pi = P{ > 0 such that 

A'nPx + P i(Aii)*<0Vz. 

After computing such a matrix Pi, define matrix W\ = (3 Pi where /? > 0 is a positive 
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scalar. Next,, we shall show that there exists a 0 such that (5.11) holds for a matrix 
11 ! = 0P i . for this purpose, define a set of positive definite matrices {QJ, % = 1, 2, . . . , l 
whereQ, = -[,l l 11 P, + P 1 (^ u )*]. 

By substituting this W { = 0P X in (5.11), we obtain 

P[-Q* + PPiPi ] < 0 V t (5.12) 

It, is evident, from the above expression (5.12) that each matrix -Q { is a negative definite 
mat rix and 1\ 1\ is a positive definite matrix. So by taking sufficiently small positive value 
of <i, (5.12) will be satisfied. The value of 0 is computed as follows. 

Foi each index t, define the scalar (3 1 such that the following inequality holds 

A[-Q, + APiPi]<o. 

One value of such /i, is computed by the formula 

A < A m [Q l ]/A M [P 1 P 1 ] 

where A m [.| is the minimum eigenvalue of [.]. The required value of 0 to satisfy (5.12) is 
selected as 

0 = min{A},* = 1,2 ,...,£. 

This completes the proof of theorem. □ 

Theorem 5.4 A act of partially normal TD systems {A t , B, A h ] is simultaneously quadrat- 
irallij stalnii zablc. 

Proof : Using the lemma 3.3 (page 40), the proof is similar to that of theorem 5.3. □ 
Remarks 5.3 : So far we have concentrated on designing state feedback controller for 

gy 

TD systems. The static output feedback controller design can be carried for TD systems 
in the same way ;is the results of chapter 3. However we have not discussed these issues in 


(let ail. 
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5.4 Controller Design for Xime Varying System 

In this section < onstant feedback controller design for time varying system is presented. 
In the previous chapters, we have identified different classes of systems for which SQ 
stabilization pioblem is solvable based on Lyapunov stability theory. Here we derive the 
stability results foi different classes of time varying systems using the SQ stabilization 
results of previous chapters. This is discussed below 
Consider a linear time varying system ( A(t),B(t )) described by 

x{t) = A(t)x(t) + B(t)u (5.13) 

The problem considered in this section is to design a static state feedback controller u = 
— 1\ x (t ) such that the equilibrium state of the closed loop time varying system 

x{t) = ( A{t ) - B(t)K)x(t) (5.14) 


is exponentially stable. 

(liven two sets of constant matrices {.4 t },i = 1,2, ...,Z and {£j,r = 1,2 define a 
parametrized family of systems described by state space equation 

i i 

(t)A t x(t) + Y Oi l (t)B l u = A a x + B a u (5.15) 

trA 1=1 


where n t (t) > 0 is a piecewise continuous scalar function and £(=i <* t (t) > 0 for all t > 0. 


0 


In this case matrix B x satisfies assumption 2.3, i.e., B x = BG t where B = 

(G\ i O'*) > 0. This parameterized family is represented by the following set S as 

t i * 

S {(.4„, B tl ) : A n = Y B a = Y a ^) B ^ ^ 0 V i and £> t (0 > °>- 


and 


.=i 


For different instant of time, nr,(0 takes different values. Suppose, a time varying system 
(.4(Q, £(0) 6 S for all time t, then matrices A{t) and B{t) are written as 


A(t) = Ya l (t)A i , B(t) = Ya t (t)BAt 

X=l 

The question is whether there exists a single control law u = -Kx such that the time 
varying systems (5.13) along with the control law is exponentially stable. To give an 
answer of this question, we next introduce the following definition. 
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Definition 5.2 : I lie linear time varying system (5.13) along with a control law 
u is exponentially stable if there exists a quadratic Lyapunov function f/(x) = 

jit, with / 0. such that the V (x) < 0 along the solution of closed loop system. 

Remarks 5.2 : For exponentially stable time varying system, the solution of closed loop 
system (5. Id) will converge to 0 as t -4 oo [73]. 

Next, wo shall state* the stabilization results for time varying systems. 

Theorem 5.5 The. time, varying system (5.13) along with a state feedback control law 
ii -. ~K.r is exponentially stable if the system realization ( Aft ), B(t)) e S V t > 0 and the 
matrices { . 1, } and { If ) satisfy any one of the following conditions. 

(a) { 4 , If } are matched systems 
or 

(b) {.dj, If] arc simultaneously transformable into H ess enberg form 
or 

(c) {.-1,, If} are partially commutative, or partially normal family. 

Proof : The theorem can be proved by employing the results of chapter 2 and chapter 3. 
Hen* we provide the proof of the part (a) and the proofs of other parts will be similar, 
proof of (a) : Since the systems {.4 t , B t } are matched systems then there exists a 
controller gain A and matrix \V > 0 such that the following inequalities hold by theorems 

2.3 and 2.7 

(.i, - ifKyw - 1 + vv '- 1 (4 - Bjcy < o v * = 1, 2 , . . . , i (5.16) 

Now consider the quadratic Lyapunov function V{x) = x*W~ l x. The derivative of this 
function along the solution of closed loop system is 

f(.r) x*(Ol(.-U^-*D(t)A')nT- 1 + ir- 1 (A(t)-R(t)AlM£) 

- x*(f)(E «,(«)[ (4 - BJQ’W- 1 + VF _1 (4 - B t K)* })x(t) 

Using (5.1(1) and noting that positive combination of negative definite matrices are negative 
definite, \'(x) < t) V t > 0. This proves the exponential stability of the time varying 

system (5.13) with a state feedback control law. □ 
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Definition 5.2 : I ho linear time varying system (5.13) along with a control law 
11 * .pnnt nt i, illy stable if there exists a quadratic Lyapunov function V{x) = 

xU>r ' wi,h P surh that thc V{x) < 0 along the solution of closed loop system. 

Remarks 5.2 : For exponentially stable time varying system, the solution of closed loop 
system (5.1 1| will converge to 0 as t -» co [73], 

Next we shall state the stabilization results for time varying systems. 

Theorem 5.5 I he turir. varying system (5.13) along with a state feedback control law 
u ~ Kx is exponentially stable if the system realization ( A(t),B(t )) € S V t > 0 and the 
matrices {.1,} and { If } satisfy any one of the following conditions. 

(a) arr matched systems 

or 

(h) {.•!,, If } are simultaneously transformable into Hessenberg form 
or 

(r) {/I,, R,} are parti ally commutative or partially normal family. 

Proof : The theorem can be proved by employing the results of chapter 2 and chapter 3. 
Here we provide the proof of the part (a) and the proofs of other parts will be similar, 
proof of (a) : Since the systems are matched systems then there exists a 

controller gain K and matrix W > 0 such that the following inequalities hold by theorems 

2.3 and 2.7 

(.4, - B t Kyw~ l + W~ l {A t - B t ky < 0 V * = 1, 2, . . . , l (5.16) 

Now consider the quadratic Lyapunov function V(x) = x*\V~ l x. The derivative of this 
function along the solution of closed loop system is 

V{x) = x'it^A^-B^KyW^ + W-^A^-B^K^xit) 

= x*(f)(X>(i)[(A - BiKYW- 1 + W~ l (A l - BiK)*])x(t) 

1=1 

Using (5.16) and noting that positive combination of negative definite matrices are negative 
definite, V(x.) < 0 Vi > 0. This proves the exponential stability of the time varying 
system (5.13) with a state feedback control law. □ 
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5.5 Conclusion 

In the first part of this chapter, the controller design problem for a collection of time delay 
system is discussed. Different classes of delay systems are shown to be stablizable by a 
single static state feedback controller. These results are valid for matched uncertain system 
as well as certain class of mismatched systems. These mismatch classes of systems can be 
thought of as a major contribution in this direction compared to the existing literature 
where the existence conditions of only matched time varying uncertain systems have been 
ensured. These results can be extended using memoryless output feedback controller using 
the results of chapter 3. 

Deviating from time delay systems, the exponential stability results of linear time vary- 
ing systems are considered in section 5.4. The SQ stabilization results of LTI system is 
successfully utilized to design a single constant feedback controller for special classes of 
time varying systems. 



Chapter 6 

Applications to Case Studies 


In this chapter, theory developed in this thesis has been applied to some practical systems. 
For this purpose, we consider two examples, namely Aircraft control and Robot trajectory 
control. In the first part of this chapter, we discuss the feedback controller design of an 
aircraft model which is uncertain. In the second part of this chapter, the trajectory control 
of robot manipulator has been discussed. 


6.1 Robust Stabilization of an Aircraft 

Consider the 4th order model [40] of the longitudinal motion of an aircraft described by 


x = A(r)x + Hu 


where A(r ) 


0.0 

1.0 

0.0 

0.0 

1.4010 x 10" 4 

r 

-1.9513 

0.0133 

-2.5050 x 10 -4 

1.0 

-1.3239 

-0.0 

-0.561 

0.0 

0.3580 

-0.0279 


H 


0.0 


0.0 


-5.3307 6.447 x 10" 3 

-0.1600 -1.155 xl0“ 2 
0.0 0.1060 


0.0 

-0.2667 

-0.2511 

0.0862 


( 6 . 1 ) 


and 


The state vector are defined as x — ^ £i x 2 x 3 £4 j control input as u 


Ui u 2 u 3 
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The state variables and control inputs are defined as follows : 
X\ : pitch attitude 
x 2 : pitch rate 
x 3 : angle of attack 
x.i : air speed 
ui : elevator angle 
u-z ' throttle control 
u-,i : flap deflection. 


In the above equation (6.1), A(r) matrix is uncertain and H matrix is fixed. The 
uncertain element r varies from —0.558 to -3.558. 

Define new set of control inputs v = Gu where the nonsingular matrix G is given by 
-5.3307 6.447 x 10~ 3 -0.2667 

-0.1600 -1.155 x 10" 2 -0.2511 

0.0 0.1060 0.0862 
Using this new set of control inputs, we define a set of two extreme plants {A, B } described 

<iS 

x = A,x + Bv (6-2) 


G 


where B = 

A 2 = -4(r) when r 2 = 
general family of systems 

33) we compute the Lyapunov matrix as W = 


Using the algorithm of section 2.9 (page 


1 

-1.1 

0.5 

0.35 

-1.1 

2.46 

-0.35 

-0.235 

0.5 

-0.35 

1.41 

0.295 

0.35 

-0.235 

0.295 

1.212 


In the above equation matrix A\ = A(r) when r — —0.558 and 

3.558. These two systems may belong to matched system or a 
defined in chapter 3. 


Using the model { A U B }, the stabilizing controller gain is computed as 


K = 


1.2729 1.0417 -0.1667 -0.1250 

-0.7067 -0.1667 1.1167 -0.1000 

-0.4988 -0.1250 -0.1000 1.1750 
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Figure 1: Distribution of eigenvalues of A{r) - HK with variation in r 


The actual control law for the stabilization of {A t ,H} is given by u — G l Kx where 


K = 


-0.4193 

-7.4912 

3.4260 


-0.2431 

-1.9167 

0.9068 


0.2727 

2.9237 

-4.7554 


0.0443 

11.203 

-0.1453 


It must be noticed from (6.1) that single uncertain element r varies between -0.558 and 
-3.558. The actual uncertain system ( A(r),H ) can be written as the convex combination 
of two systems (A\,H) and ( Ai,R ). More specifically, 


2 1 

A(r) = ot % A u H = OiiH, a t > 0, ]T) CKi = 1. 

»=i 1=1 

So the uncertain system (A(r),H) is stabilizable by a single control law u = -Kx. To 
verify this, we plot in figure 1 the distribution of the eigenvalues of matrix (A(r) - HR) 
with the variation of r in the range -0.558 to -3.558. This confirms that all eigenvalues 
of matrix (A(r) - HK), with -3.558 < r < -0.558, have negative real part. 
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6.2 Controller Design For Robot Manipulator 

It, has boon shown in chapter 2 that a set of systems in block companion form is simul- 
taneously quadraticallv stabilizable by a single static state feedback controller. Based on 
this result, a new approach to robot tracking controller design is presented. The proposed 
control scheme consists of a feedforward controller based on the inverse dynamics of the 
robot and a feedback controller. The non-linear model of the robot is viewed as piecewise 
LTI systems obtained by linearizing the model at selected number of points on a specified 
trajectory in the joint space. The collection of all the LTI systems constitutes a set in 
which each member is observed to be in block companion form. For such class of svs- 
tems, an algorithm for the design of a single stabilizing feedback controller is presented. A 
numerical example of a two link manipulator is considered to validate the proposed theorv. 


6.3 Preliminaries 

Present-day control literature reveals umpteen number of research works on the tracking 
control of non-linear systems. Of the many approaches available, the application of linear 
multi-variable theory has received wide acceptance because of its simplicity and ease of 
understanding. In this approach, the non-linear system is viewed as a set of LTI systems 
over some specified trajectory. In such a situation, it is always desirable to have a single 
stabilizing controller for the above mentioned set of LTI systems and hence for the overall 
non-linear system. This motivates the application of the simultaneous stabilization problem 
which amounts to the design of a single feedback controller for stabilizing a finite number 
of plants. 

The problem of control of robot manipulators for trajectory tracking has long attracted 
the attention of researchers and engineers (see, for example [37, 63, 19] and the references 
therein). This problem is concerned with the design of an appropriate control scheme to 
track time varying joint space or task space trajectories. In this chapter we shall focus 
only on the former. There are, at present, numerous approaches to the design of robot 
controllers of which the Computed Torque Technique (CTT) is most common. The CTT 
scheme is based on the global linearization of the robot dynamics by means of a non-linear 
feedback controller. However, the order of complexity of the controller is same as that 
of the robot model [52] which makes it computationally expensive. Seraji [60] proposed 
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a simple linear control scheme using a feedforward and state variable feedback control, 
llie robot dynamics is approximated by a set of linearized models over the trajectory. 
However, the control scheme requires the updating of feedback gains which can be realized 
by off-line or on-line computation with gain scheduling. This process of updating results 
in the overall system to be time varying and may lead to deterioration of performance 
duo to switching of the gains depending on its rate. Moreover, for off-line gain scheduling, 
there is a trade-off between the performance of the controller and the memory requirement 
for storing the the pro-computed gains. Similar control scheme has also been proposed by 
Swamp and Gopal [66] using a performance optimization with eigen structure assignment. 
They also considered the uncertain dynamics of a manipulator to solve tracking problem. 

In chapter 2. we have presented an algorithm for the construction of a single static 
state feedback controller for a set of LTI systems This feedback controller is designed 
based on the construction of a single Lyapunov function for the set of LTI systems and its 
applicability to time varying systems is discussed in chapter 5 (section 5.4). The theory 
developed so far is then used to design a feedforward and a single static state-variable 
feedback controller for a serial chain robot manipulator. The feedforward controller ensures 
the tracking of the manipulator joint trajectory while, the feedback controller ensures the 
exponential stability of the overall closed loop dynamics. The proposed approach has a 
distinct advantage over the previous approaches [60] in that we do not require any feedback 
gain scheduling. 


6.4 Application to Robot Control 

In this section, the theory developed thus far are applied to the problem of robot control. 
The problem considered here is to develop a control scheme which ensures that the joint 
angle vector 9(t) tracks any desired (reference) trajectory 0 r (t) where 9 r (t) is a vector of 
smooth time functions. A simple linear control design consisting of a feedforward and a 
static state- variable feedback control is proposed. 
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6.4 Application to Robot Control 


6.4.1 Robot Dynamics : Linearized Model 

Using the Lagiunge-Eulei formulation, the equations of motion of a serial-chain manipu- 
lator can he derived in the following form [19, 37, 63] 

M{e)e + N(e,e) + G(e) = T (6.3) 

when; M{9) is n x n symmetric positive definite inertia matrix, N{9,9) is n x 1 Coriolis 
and centripetal forces, G(0) is n x 1 gravity loading vector and T is the n x 1 vector 
of joint torques. 0,9,6 are n x 1 vectors of joint angles, velocities and accelerations, 
respectively. The elements of M(9), N(9,6) and G{9) are highly non-linear functions of 
the joint configurations, velocities and the payload. 

After tin; desired trajectory 9 r (t) has been determined in the planning stage, the ap- 
propriate torque is required to result in the desired motion. Under the ideal conditions 
(the model is accurate and the absence of any disturbances), the desired torque is given 
by [37] 

T d = M(6 r )9 r + N{6 r , 0 r ) + G{9 r ). 

If this torque is applied to the manipulator modeled by (6.3), the resulting motion tracks 
the desired trajectory under ideal circumstances. Because of disturbances and inaccuracies 
in the system model, the actual motion of the manipulator usually deviates from the desired 
trajectory. These undesirable deviations can be corrected by means of feedback controller. 
So using PD type feedback controller, the total torque which can be used for the tracking 
purpose is 

T = T d + K p (9 r -9) + K v (9 r -9). 

where K p and I\ v are the appropriate feedback gain. 

Since the actual motion of manipulator deviates from the nominal trajectory, it is 
necessary to study deviation versus time behaviour from the nominal trajectory and the 
compensation for perturbations. These goals can be achieved by using linearized model 
of the manipulator dynamics that describes the small perturbations about the nominal 
trajectory. The determination of such linear model is then discussed which will be used 
for the purpose of controller design. 

In the following developments, a specified trajectory in the joint space is considered as 
a collection of points lying on it. Consider the i‘ h nominal point (node point) denoted by 
S< = The nominal joint angle vector 0 m and its higher derivatives are 
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taken on the prespeoified trajectory and T m is the corresponding nominal torque. At the 
time of linearization, it has been assumed that the inertia matrix does not change appre- 
ciably for a small deviation x about the nominal point, i.e., M(9 ni + x(t)) « M(6 m ) [60]. 
The linearized dynamic model, with respect to the i th nominal point, may be represented 
in the following form 

N 2l x + N u x + N 0l x = u (6.4) 


where N 2l , N u , Ma an; n x n constant matrices obtained as 


N 2i = {M(9)} s >,N u = 


8N' 

86 


, N 0l = 


d(N + Gj 

89 


The vectors x and u are deviations from the nominal values of joint positions and torques, 
respectively. Defining the state variables as (x,x), the above linearized equation (6.4) is 
represented in the state space form as 


d_ 

dt 


x 

x 


0 In 

-n£n 0 . 


-] 

X 


0 



+ 



X 


1 

-i 

£ 
i 


u. 


(6.5) 


It may be noted that the robot model over the trajectory is a slowly time varying system [19. 
CO] and is approximated by piecewise (in time) LTI in nature. The matrices N 2i ,N u ,N 0t 
are evaluated at different node points * = 1,2,... and are used to obtain a set of LTI 

systems in the form of (6.5). 


6.4.2 Controller Design 

In this section, a control scheme is presented based on the linearised model obtained m 
(6.4). As may be noted from (6.5), the linearized dynamics is always controllable and the 
matrix jVp > 0. ,V 2 , being the inertia matrix of the robot evaluated at a nominal point. 

At a nominal point on the trajectory, all the non-linearities in the robot dynamics 
are exactly canceled by the nominal joint torque T„, since it is calculated based on the 
non-linear inverse dynamics (6.3) at the nominal point. About this point, the linearized 
dynamics is represented by (6.4). The proposed control scheme consists of feedforward and 
feedback controllers based on this linearized model. The implementation of the scheme for 
the non-linear dynamics of the robot is shown in Fig. 2 where s is the differential operator. 
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Figure 2: Structure of Robot Control Systems 

The total control law is given by 

U ( t ) = u f f(t) + Ufb{t) + T m = u(t ) 4- T ni (6-6) 

where. and ujb(t) are the feedforward and feedback controllers respectively and 

U{t) = Ufj(t)+Uf b {t). 

The feedforward controller, computed based on the inverse dynamics of the linearized robot 
model (6.4), provides tracking ability and is given by [60] 

Uff(t) = ^2iX r + NfoXr + NfoXr. ( 6 . 1 ) 

where x r denotes the deviation of the nominal point on the reference trajectory, i.e., the 
reference vector 6 r (t) = 6 ni +x r . 

The feedback controller provides the asymptotic stability of the overall system and is 
discussed below. Substitution of the expressions for u(t) in (6.4) yields the closed loop 
dynamics as 

N 2 & + NuX + N&X = N2iX r + Ni t X r + NoiX r + Ufb(t) 

=> iV 2l e + N\ x e + iVo,e = ~ u fb(t) 
where. e(t) — x T (t) — x(t) is the tracking error. 


( 6 . 8 ) 
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Introducing the' state' vector r) = [ e(tY 


e(t)* ]*, (6.8) can be written as 


v{t) — + B t Ufb(t ) 


(6.9) 


where?, ,4, 


--A^/Vo, -N£N lt 


and B, = 


0 

~ N 2i l 


As the robot, manipulator moves over the reference trajectory, the A z and B l matrices 
e-liange' continuously resulting in a linear, slowly time varying system [19]. Hence, on a 
discretized trajectory as in the present case, the linearized dynamics may be approximated 


b\ piee ewise constant (in time) A x and B l matrices [60]. Asymptotic tracking is achieved 
here; in the sense that r/(£ ) -4 0 as t -4 oo. 

At this point a discussion on the previous research is in order. In the existing ap- 


proaches, the state feedback controller of the form Ufb(t) — K x r](t) has been considered. 
Ihe matrix A , is computed at every nominal point S 1 on the reference trajectory such 


that, all eigenvalues of the closed loop matrix {A, + B l K l ) are in the open left half plane 
and remain unchanged. The gain matrix K t is updated successively as the manipulator 
moves. This calls for a scheduling of the feedback gains which may be computation in- 
tensive and may deteriorate tracking performance depending on the switching frequency 
[60, 70]. This motivates us to propose a new design methodology in which a single feed- 
back controller is designed to ensure the exponential stability of the tracking error dy- 
namics of the robot over a specified trajectory. Using a state feedback controller of the 
form Ufb(t) = Kr](t) = K p e(t ) + I\ v e(t), the overall closed loop error dynamics over the 
trajectory is of the form 


fl(t) = AjXyqit) 


( 6 . 10 ) 


where, A c {t) = (.4, + B l K),t t < t < t l+] _. The time instants t x and t l+i correspond to 
two consecutive nominal points. 5, and S t+ i, respectively. The key steps in the design 
procedure are : 


• Finding a feedback controller of the form u/&(£) = Kr]{t). 

• Finding a quadratic Lyapunov function of the form v(rj) = rj*Pr] such that v(rj) < 0 
along the solutions of (6.10). 

It is observed from (6.9) that the error dynamics obtained from the linearized model (6.5) 
is in block companion form. The collection of all such systems, defined on the trajectory, 
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constitutes a sot of LTI systems {A u B t }. It is also noted from the expression of B t matrix 
of (6.0) that the matrix > 0, N 2l being the inertia matrix of the robot evaluated 
at a nominal point. I his allows us to use the result of theorem 2.5 and the subsequent 
algorithm for designing a single state feedback controller gain K. Further, a common 
Lyapunov matrix P > 0 exists such that 


(-4, + B l K)*P + P(A t + B X K) < 0, V i. (6-11) 

Next, to show the exponential stability of (6.10), consider the Lyapunov function v(rj) = 
if Pi). Taking the time derivative of v(rj) and using (6.11) yields 

v(rj) = r)*[A c (t)*P + PA c (t)]r) < 0 Vi. 

Thus, the asymptotic tracking is achieved in the sense that the tracking error 77(f) -4 0 as 

t ->■ oo. 


Remarks 6.1 : In the above discussions, the exponential stability of linearized error 
dynamics is proved. However, the error dynamics for actual nonlinear systems is also 
stable locally with respect to the given trajectory. 

Remarks 6.2 : In this example, we have not considered the uncertainty of parameters or 
load. The complete investigation should be done in this direction. However, our motive 
behind taking this example is to show the application of SQ stabilization results to a class 
of linear time varying systems. 

6.4.3 Numerical Example 

In order to validate the theory presented, a numerical example of a two-link planar ma- 
nipulator, as shown in Fig. 2, is considered. The inertial and geometric parameters of 
the manipulator are taken from [60]. For the purpose of simulation, the manipulator is 
assumed to move from an initial joint angle configuration of [tt/ 2, 0] to the final configu- 
ration of [0, tt/ 2]. The initial and final joint angle velocities are assumed to be zero. This 
is achieved by planning the reference trajectory in the following functional form 

d rl {t) = 7r/2[3ezp(— f/0.3) — 4ezp(— f/0.4)], 

d r2 (t) = ir/2[l + 3exp{-t/0.3)-4exp(-t/0A)\ 
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where t varies from 0-5 sec. The exponential form of the trajectories are chosen to ensure 
smoothness. The reference trajectory is discretized by choosing 50 equidistant node points 
at which the dynamics is linearized to obtain a set of LTI systems. Following the algorithm 
presented in section (2.9), a single state feedback controller is designed as 


-87 43 10.43 -124.0 -0.29 

-36.94 -25.52 -0.293 -123.5 


The feedforward control law are calculated on-line using (6.7) and updated at each node 
point on the reference trajectory The simulation was performed by integrating the non- 
linear equations of motion of the manipulator (6.3) using the fourth order Runge-Kutta 
algorithm. The reference and the actual joint angle trajectories are presented in Figs. 3 
and 4. The small joint angle errors over the trajectory shown in Fig. 5 indicate the efficacy 
of the proposed control scheme. 


6.5 Conclusion 

In this chapter, we have applied the results of SQ stabilization problem to design of a 
feedback controller for an aircraft model and a Robot model. In the first case, we consider 
an uncertain aircraft model to design a single state feedback controller. In the second case 
we consider the tracking problem of a robot manipulator. In the latter case, the feedback 
design problem is thought as a design of a feedback controller to stabilize a time vary- 
ing svstem. The key feature of the proposed approach is that it obviates feedback gciin 
scheduling. This feature is advantageous from both the performance and the implementa- 
tion points of view. 








Chapter 7 
Conclusions 


7.1 Conclusions 


An important problem in control theory is to stabilize finite number of linear time invari- 
ant systems by a single controller. To fulfill this requirement, this thesis considers the 
stabilization of a finite set of systems by a constant feedback controller. The controller is 
also ensured a common quadratic Lyapunov function for all the systems. The stabiliza- 
tion problem of a set of LTI systems is formulated as a joint search problem for a single 
feedback gain and a quadratic Lyapunov function and referred as simultaneous quadratic 
stabilization problem. A number of parametric conditions describing the set of linear sys- 
tems in state space are determined. These conditions guarantee the existence of constant 
feedback controller for which simultaneous quadratic stabilization problem is solvable. For 
this purpose, we consider a set of LTI systems {A t , B t } described by state space equation 
where 


A = 




A 1 


21 


12 
A’,, 


,B X = BG t , B = 


0 

Im 


and ( Gi + G*) > 0. The classes of systems 


which are simultaneously quadratically stabilizable by static state feedback controller are 
as follows : 


• Matrices A‘ u and Aj 2 are same for each system. 


• Each system matrix A x are in Hessenberg form. 


• Matrices {A* u } are Hurwitz and commutative pairwise. 

• Matrices {AjJ are Hurwitz and normal. 
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• Matrices {A\ 2 } are quasi sign invariant. In this case the number of inputs m are 
constraints as follows : 

— m > nj 2 if n is even and 

— m > (n + l)/2 if n is odd. 

» 

The above classes of systems (except first one) are different from the matched uncertain 
systems. It has also been shown that if a set of system is stabilizable by a state feedback 
then their positive combination will also be stabilizable. 

To design static output feedback controller, we consider a set of square (equal number of 
inputs and outputs) minimum phase systems {4, B u C}. Under the assumptions that the 
high frequency gain is non zero, the following conditions are derived for the SQ stabilization 
of (4, Bt, C} by static output feedback controller : 

• The matrices A\ x and A\ 2 are same for each system. 

• The matrices (4i ~ A\ 2 C 2 l C\] are Hurwitz and pairwise commutative. 

• The matrix (A l u — A\ 2 C 2 l C\} are Hurwitz and normal. 

The SQ stabilization problem is also extended for special class of time delay systems 
{A t ,B,Ah}. The existence of a stabilizing controller is assured for this class of systems. 
The SQ stabilization problem is also shown to be useful for the stabilization of time varying 
systems. As an application, we design controllers for an aircraft and a robot manipulator. 

In this thesis, different class of systems in state space framework are identified for 
which SQ stabilization problem is solvable. These results are in fact applicable to rather 
specialized situations. However, it is true that if a finite set of systems is SQ stabilizable 
then the set of systems, consisting the neighbourhood of each system, is also SQ stabilizable. 
In other words, if (4, #i) and (Ao, Bo) are SQ stabilizable then there exists neighbourhood 
of each plant Ni and N 2 such that the systems belonging to AfiUA^ are also SQ stabilizable. 
Based on this observation, the results of this thesis will be useful for the robustness analysis 
of general class of systems. 
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7.2 Scope of Further Research 

In this thesis, we characterize different classes of systems for which SQ stabilization problem 
is solvable Once the existence conditions are satisfied, the controller gain is computed from 
the algorithm. However, it is not ascertained that whether controller gain K is of minimum 
norm. So a procedure has to be found out to compute SQ stabilizing controller gain of 
minimum norm. This will be useful to avoid the saturation of sensors. 

In practice controller should ensure certain other performance criteria in addition to 
stability So it will be interesting to design a SQ stabilizing controller which will ensure 
certain performance robustness. 

In recent times, the ideas of switching between different controllers are used in intelligent 
control and applied to complex systems in the presence of structural failures [48, 46, 47]. In 
this approach, it is assumed that the system can be in only a finite number of configurations 
after failure. Each configuration corresponds to a compact set of points in the space of 
system parameters. One way of reducing the number of switching is to design a single 
controller for each configuration. Thus after determining a configuration to which the plant 
belongs (using pattern recognition methods) , it is necessary to design a single controller for 
each configuration. At this stage, SQ stabilization will be useful to determine controllers 
as well as analyze the stability of such switching systems. 



Appendix A 
Long Proof 


A.l Proof of Lemma 2.3 


Proof: From the results of controllability [79], if (A t , B ) is controllable then for any matrix 
H , (-4! + BH, B) is also controllable. Using the partitioned matrices of (A, B) given in 
(2.12), the following simplifications are made. 


( A n A\ 2 

\ . Ah A. l 22 


L 


is controllable. 




4* 

•“M2 

0 


0 

4n 

0 


[A 


■21 /1 22b 


0 


is controllable. 


is controllable. 


=> Rank 

=> Rank 


0 A\ 2 A\ x A\ 2 {A l n ) 2 A\2 ... _ 

Im 0 0 0 ...J ~ 

A\ 2 A' n A\ 2 {A\ x ) 2 A\ 2 ...] = n — m. 


=» (A l u ,.4' 12 ) is controllable. 

Similarly, we can show that if a system (.4,, B) is stabilizable then the system (A 21 , A x2 ) is 
also stabilizable. This can be proved using the eigenvalue eigenvector criteria to test the 
stabilizability [2, 34] and is discussed below. 
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The stabilizability of the system (A,, B) is equivalent to the following [2, 34] : 
x* [A I n — A{ B] = 0 for some x* ^ 0 implies FeA < 0. 

Since ( A t ,B ) is stabilizable, so 

x* [A I n — A z ] = 0 and x*B = 0 for some x* ^ 0 implies Re A < 0 . 

Note that here 0 is a vector or matrix of appropriate dimension. The vector x which sat- 
isfies x*B = 0 is given by the form 



where z £ R n m . Using this x , the following simplification is made 


x* [A/ n — A] = 0 
=> Z* [A In-m ~ A l n A\ 2 = 0 

Since (A t , B) is stabilizable, so Re A < 0. Hence it is evident from the above that (A J n , A[ 2 ) 
is also stabilizable. □ 


A. 2 Proof of Theorem 2.5 

Sufficient: Suppose there exists a matrix W\ = W* > 0 and a matrix F such that the 
following inequalities holds, 

[(All + A\ 2 F)W x + WMn + ^12^)1 < 0 V 2. 

We need to find a positive definite matrix W and a scalar 7 > 0 such that $*, i = 1 , 2, . . . , l 
will be negative definite where Tj is defined by 

= [A t W + WA; - (7/2 )B{Gi + G:)B*}. (Al.l) 

For this purpose we construct a positive definite matrix W in the following way. By 
using the matrices F and W\, choose W 2 = FW'i. Then one can select W3 such that 

[" \v, W 2 1 

W 3 - W:Wy X W 2 > 0 which ensures that W = is positive definite matrLx. 

\ V 2 v* 3 

Such a choice of W 3 is always possible. 

With the above matrices W and F, the matrix is partitioned as 

$ t = [ ^ ^ with <f> l n e R*-™* n - m . 

(<^12)* 4 > 22 . 
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The different submatrices are defined as follows (by noting the special structure of ma- 
trix B defined in assumption 2.3, page 32) 


4>u = (X* u + A‘ u F)W l + WMn + A\ 2 F)' 

<t >\2 = A' u w 2 + a\ 2 w 3 + w^y + w 2 (A‘ 22 y 
= A‘ 21 w 2 + A' 22 w 3 + + w 3 (Ai 2 y - 7(G, + G’) 

In the above equation F = W^W^ 1 . 

The matrix will be negative definite if [38] 

(/) 0* n < 0 and (II) 022 - (0i2)* (0ii) _1 0k < °- 

It is observed that the matrix 0 z n is negative definite for all i by the hypotheses of 
theorem. Also it should be noted that 0* u and 0k are independent of 7 and the matrix 
(G l + G*) > 0. So the matrix 0 22 can be made negative definite by choosing a sufficiently 
large value of 7. Since d> l n and 0* 2 are independent of 7, we can choose 7 such that 

022 - (0k)* (0ii) _1 0k < 0 Vi. (A1.5) 

The value of 7 is computed as follows. Using the expression of 0k from (Al.4), the above 
inequalities (Al.5) are simplified as 

lIm(Gi + <?*) > S l — (0k)* (01 1 ) 1 012 ^ * (A1.6) 

where S, = Ak^k + AkUk + ^2(^21)* + ^(Ak)*- 
Now, for each index 2, choose a scalar 7, such that 

7i > -W[ (Si — (0k)* (0ii) 1 0k) (^ A G t ) ] 

where A,u[.j is the maximum eigenvalue of matrix [.]. 

Let the value of 7 = max^J, 2 = 1, 2 , . . . , L For this value of 7 and \\ , the matrix 
defined in (.41.1), will be negative definite for each 1 = 1, 2 , . . . ,1. The controller gam 

then will be K = ( , y/2)B*W~ l . 


(A1.2) 

(A1.3) 

(Al.4) 
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Necessary : Suppose, there exists a scalar 7 > 0 and a matrix W > 0 such that the 
following matrix inequalities hold : 


= [A t W + WA; - {y/2)B(G t + G*)B*} < 0 Vi. 


(A1.7) 


Substituting B = 

tioned as 
= 


0 


and partitioning W = 


Wi W 2 
W 2 * W 3 


, the matrix are parti- 


. By using F = W 2 W X 1 , the different submatrices are defined as 


4>u 4>n 

( 012 )* *^22 

in (Al 2-A1.4). Since the matrix is negative definite for each 1 , so matrix <^ u is also 
negative definite for each 1 . Hence the proof. □. 
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Long Proof 


B.l Proof of Theorem 5.1 

Suppose there exists a matrix W x — W{ > 0 and a matrix F such that the following 
inequalities hold, 

l(A\ x + A\ 2 F)W x + W X {A l u + A\ 2 F)* + W x {I n . m + F*F)W X ) < 0 V * 

We need to find a positive definite matrix W and a scalar 7 > 0 such that 1 — 1, 2, . . . , l 
will be negative definite where <E> t is defined by 

$ t = [A % W + WA* - iBB* + A h A* h + WW]. (B2.1) 

For this purpose we construct a positive definite matrix W in the following way. By 
using the matrix F and W x , choose W 2 = FW X . Then one can select W 3 such that 

I" W 2 ] 

Wo - WJWrHV-, > 0 which ensures that W = is positive definite matrix. 

2 1 “ [ ir 2 * W 3 _ 

Such a choice of W 3 is always possible. 

With the above matrix W and F, the matrix is partitioned as 

q <t>\ 1 <Pn with p € R n-mxn-m 

[W <P22\ 
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The different submatrices are defined as follows by taking care of special structure of 
B and Ah matrices given (5.9) (page 64) 


fix = {*U + ^12 F)Wx + + A\ 2 FY + Wifo-n + F*F)Wi (B2.2) 

<t>\ 2 = a\ x w 2 + a\ 2 w 3 + wm i 2 iY + W 2 (A i 22 y + [w x w 2 ][w; w 3 ]* (B2.3) 

(f >*22 = A' 2X W 2 + A 22 \V 3 + W*(A l 2i y + W z (A' 22 y - 7 I m + DD* + [W 2 * W 3 ][W^ W 3 ]\ (B2.4) 
The matrix <]?, will be negative definite if [38] 

(/) 0- u < 0 and (II) - yy («,)-■ < 0. 

It is observed that the matrix <f> xl is negative definite for all i by the hypotheses of theorem. 
Also it is noted that <j>\ x and <j>\ 2 are independent of 7 and the matrix (G t + G*) > 0. So 
the matrix (j> 22 can be made negative definite by choosing a sufficiently large value of q. 
Since 4> l n and <j>\ 2 are independent of 7, we can choose 7 such that 

<*22 - Ma)* Wn)' 1 fv 2 < 0 V z. (B2.5) 

The value of 7 is computed as follows. Using the expression of cj) 22 from (B2.4), the above 
inequalities (B2.5) are simplified as 

1 > — (<^12) (0ii) 1 ^12 ^ ® (B2.6) 

where S* = (A\ X W 2 + A^Wz + W;(A' 2X )* + W 3 {A\ 2 )* + DD* + [W 2 * W 3 ][W . 2 * W 3 ]*. 

For each index i, choose a scalar 7 , such that 

7i > -W( Si ~ (^12) (01l) 1 012 ) 

Let the value of 7 = max { 7x } , 2 = 1,2, For this value of 7 and W, the matrix 

$ t , defined in (£2.1), will be negative definite for all i. The controller gain then will be 

K = (7/2)£*ir- 1 . □ 
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Quadratic Form of a Matrix 

In this appendix, a few results regarding the quadratic form of a matrix are presented. 
However, details are available in standard matrix theory book (for example [30]). 

An n x n matrix M is said to be positive definite matrix iff the following quadratic 
form 

x*Mx > 0 

for all non zero vector x. 

An n x n matrix M is said to be negative definite matrix if -M is a positive definite 
matrix. 

Next the following facts are true for a positive definite matrix. 

Fact 1 : A symmetric matrix A is positive definite iff the following are true 

(a) All the eigenvalues of matrix M are positive. 

(b) All leading principal minors are positive. 

Fact 2 : Suppose that a symmetric matrix M is partitioned as M = 

and C are square matrices. The matrix M is positive definite iff 
(i) A > 0 (ii) C - B*A~ 1 B > 0. 

Fact 3 : Suppose A and B are symmetric matrices. Then the matrix A+B will be positive 

definite if 

A m [A] + Am[J5] > 0. 

Fact 4 : Suppose A and B are two symmetric and A is positive definite along with B 
positive semidefinite, then {A + B) is also positive definite. 


A B 
B* C 


where A 
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ERRATA 


The following two references will be added at the end of bibliography at the page 83. 

[84] Boyd S.P., El Ghaoui L., Feron E. and V. Balakrishnan: Linear Matrix Inequalities 
Systems and Control Theory, SIAM studies in applied mathematics, 1994. 

[85] Skelton and T. Iwasaki : Increased roles of Linear Algebra in Control Education, IEl 
Control Systems, 15(4), pp. 76-90, 1995. 

These references will be read in page 6 after the line number 7 as follows: 

Other numerical methods for solving linear matrix inequalities can be found in [84, 85 a i 
references therein]. 

The following lines are modified in the following ways. 

In page 4, line 2 from top: control system design [29]. 

In page 4, line 15 from top: reported in [53,80,42,4,35,54,68,51]. 

In page 7, line 11 from top: related results [36,9,8,12,15,74,78].... 

In page 48, line 3 from bottom: algorithmic point of view [43,20] — 

In page 71, line 9 from top: 4th order model [40,50].... 

In page 87, line 4 from top: of controllability [79,14,22,25,45,13,3], if... 



